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extension  of  the  model  to  a three-dimensional  form  for  analysis  of 
distorted  inflows,  "nie  models  satisfy  mass,  momentum  and  energy 
equations  on  a time-dependent  basis.  Compressor  stage  force  and 
shaft  work  were  determined  from  empirical  stage  characteristics 
with  corrections  made  for  unsteady  cascade  airfoil  aerodynamics. 
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CHAPTER  I 
INTRODUCTION 

Stable  aerodynamic  operation  of  the  compression 
system  of  an  aircraft  gas  turbine  engine  is  essential  for 
the  engine  to  operate  safely  and  to  satisfactorily  deliver 
the  performance  which  it  is  designed  to  provide.  The  com- 
pression system  of  an  aircraft  gas  turbine  engine  consists 
of  one  or  more  compressors  arranged  in  configurations  such 
as  those  illustrated  in  Fig.  1 (Appendix  A).^  Aerody- 
neunically  stable  operation  refers  to  the  condition  in  which 
the  compressor  delivers  a desired  quantity  of  airflow  on  a 
continuous  basis,  compressed  to  a desired  pressure  level, 
free  of  excessively  large  eunplitude  fluctuations  in  the  flow 
properties  throughout  the  compressor.  The  aerodynamic 
stability  limit  is  most  commonly  represented  as  a locus  of 
corrected  rotor  speed  points  denoting  the  minimum  corrected 
flow  rate  value  for  stable  operation  on  a plot  of  the  total 
pressure  ratio  across  the  compressor  versus  corrected  air- 
flow rate  as  shown  in  Fig.  2.  Operation  to  the  right  of  the 
curve  is  stable  and  to  the  left  of  the  curve  is  unstable. 

The  stability  limit  curve  is  often  referred  to  as 
the  surge  line  or  stall  line  of  the  compressor  stemming  from 
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the  fact  that  surge  is  one  possible  specific  form  of 
instability  experienced  in  compressors,  and  all  compressor 
instabilities  are  initiated  by  aerodynamic  stall  of  some 
blades  or  vanes  of  the  compressor. 

The  stability  margin  of  the  compressor  is  commonly 
taken  to  be  the  percentage  difference  between  the  compressor 
total  pressure  ratio  at  the  stability  limit  and  that  at  the 
operating  line  holding  corrected  airflow  constant  (Fig.  2) . 
The  position  of  the  operating  line  is  primarily  determined 
by  the  load  imposed  on  the  compressor  by  downstream  engine 
components.  However,  under  actual  flight  conditions  in  an 
aircraft,  the  positions  of  both  the  stability  limit  line  and 
the  operating  line  are  influenced  by  external  disturbances 
in  the  flow  at  either  the  entry  or  exit  of  the  compressor. 

The  overall  objective  of  this  study  is  to  develop 
an  improved  analytical  method  which  will  allow  prediction  of 
the  loss  of  stability  margin  of  the  compressor  caused  by 
external  disturbances  imposed  upon  it  by  flight  conditions 
in  an  aircraft.  Also,  the  development  of  the  method  is 
intended  to  provide  some  degree  of  physical  understanding  of 
the  processes  which  lead  to  compressor  instability. 

Prior  to  beginning  the  analysis  it  is  helpful  to 
review  the  ramifications  of  stability  on  engine  operation 
and  to  review  some  of  the  experimentally  observed  influences 
of  the  operating  conditions  and  external  disturbances  on 
compressor  stability.  Also,  a review  of  some  of  the 
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previously  accomplished  works  directed  at  stability  pre- 
dictions both  with  and  without  external  disturbances  will  be 
made. 


I.  RAMIFICATIONS  OF  STABILITY  TO  ENGINE  OPERATION 

Operation  of  the  compressor  to  the  left  of  the 

stability  limit  (Fig.  2)  will  result  in  loss  of  thrust, 

possible  loss  of  engine  control,  or  possible  engine 

structural  damage  caused  by  overheating  and  high  cyclic 

stress  induced  in  the  engine  by  unstable  flow  conditions. 

2 

Reference  [1]  presents  evidence  of  thermal  and  cyclic 
stress  effects  imposed  by  compressor  instability  for  one 
particular  engine.  Montgomery  [2]  has  compiled  a history  of 
difficulties  in  the  operation  of  several  aircraft  gas 
turbine  engines  which  have  been  experienced  because  of 
compressor  stability  problems.  Montgomery's  work  well 
illustrates  the  need  for  a method  of  determining  the  sta- 
bility limit  of  the  gas  turbine  engine  compressor, 
especially  as  they  are  affected  by  the  conditions  of  flight. 

II.  MAJOR  OPERATING  CONDITIONS  AND  EXTERNAL 
DISTURBANCES  INFLUENCING  STABILITY 

The  position  of  the  stability  limit  line  and  the 
operating  line  of  the  compressor  of  a given  design  are 


Numbers  in  brackets  refer  to  similarly  numbered 
references  in  the  Bibliography. 
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affected  in  flight  by  the  condition  of  the  flow  entering  the 
compressor  and  by  the  throttling  or  back  pressure  provided 
by  the  compressor  or  engine  component  downstream  of  the 
compressor  under  consideration.  The  entry  total  tempera- 
ture, the  Reynolds  number  at  which  the  compressor  must 
operate,  and  the  nonuniformity  and  unsteadiness  of  the  flow 
at  the  compressor  inlet  have  major  influences  on  the  sta- 
bility of  the  compressor. 

Because  of  the  variation  of  total  temperature  at 
the  compressor  entry  which  occurs  with  flight  Mach  number 
and  altitude,  the  compressor's  operating  region  moves  far 
from  the  design  point  of  the  compressor;  thus,  a relatively 
large  region  of  the  compressor  map  must  be  considered  in  any 
analysis  of  stability.  Figure  3a  depicts  the  variation  of 
operating  region  with  compressor  inlet  total  temperature  for 
a simple  single  spool  turbojet  compressor.  Figure  3b  shows 
the  influence  of  flight  Mach  number  and  altitude  on  engine 
entry  total  temperature.  The  changes  in  flight  conditions 
cause  the  Reynolds  number  to  change.  Figure  3c  gives  an 
example  of  the  Reynolds  number  influence  on  the  stability 
limit.  The  Reynolds  Number  Index  (REI)  used  in  Figure  3c  is 
the  ratio  of  the  Reynolds  number  in  flight  to  that  which 
would  exist  at  sea-level-static  conditions  at  the  same 
compressor  inlet  axial  Mach  number.  Figure  3d  shows  the 
variation  of  REI  with  flight  condition. 

Figure  4 shows  experimentally  determined  losses  in 
stability  margin  caused  by  nonuniform  (distorted)  flow 
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conditions  and  unsteady  flow  conditions  at  the  compressor 
entry  plane.  In  Fig.  4 the  loss  is  illustrated  by  downward 
shifts  of  the  stability  limit  toward  the  operating  line.  In 
some  cases,  the  disturbances  also  cause  stability  margin 
loss  by  moving  the  operating  line  upward.  For  purposes  of 
illustration,  the  entire  loss  is  shown  in  Fig.  4 as  a down- 
ward movement  of  the  stability  limit.  Figure  4a  from  [4]  is 
for  the  case  of  steady  but  distorted  flow  into  the  com- 
pressor created  by  a wire  screen  of  nonuniform  porosity 
placed  in  front  of  the  compressor.  Figure  4b  from  [5] 
illustrates  the  loss  in  stability  margin  caused  by  uniform 
but  sinusoidally  time-varying  flow  into  the  compressor 
produced  by  a rotating  flow  interrupter  upstream  of  the 
compressor.  The  case  of  combined  unsteady  and  distorted 
flow,  sometimes  called  time-variant  distortion,  is  shown  in 
Figs.  4c  and  4d  for  a turbojet  engine  compressor  [6]  and  a 
turbofan  engine  fan  [7] . The  engine  entry  flow  fields  were 
produced  by  a venturi -center body  arrangement  which  produced 
an  unstable  shock  wave  boundary  layer  interaction  which 
causes  the  adverse  flow  field  at  the  engine  entry  plane  [6] . 
Figure  43,  based  on  data  from  [8],  indicates  the  influence 
on  the  downstream  compressor  stability  limit  line  of  the 
nonuniform  and  unsteady  flow  field  exiting  from  a forward 
compressor  of  a turbofan  engine. 

In  addition  to  disturbances  in  the  pressure  enter- 
ing the  compressor,  total  temperature  disturbances  may  occur. 
Information  from  19]  shown  in  Fig.  4f,  indicates  the 
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experimentally  determined  influence  of  time-steady  spatial 
variation  of  total  temperature  entering  the  compressor. 

Rapid  uniform  total  temperature  increases  at  the  engine 
inlet  which  could  be  caused  by  ingestion  of  hot  gases  pro- 
duced by  gun  or  rocket  fire  can  also  cause  a loss  of 
stability  [10] . 

Disturbances  to  the  compressor  may  originate  at 
the  exit  of  the  compressor  as  well  as  at  the  inlet.  An 
experiment  conducted  on  a turbofan  engine  described  in  [11] 
demonstrated  the  influence  of  an  oscillating  back  pressure 
on  fan  stability. 

The  distortion  examples  of  Fig.  4 all  use  total 
pressure  or  total  temperature  variation  to  describe  the 
distortion  pattern  characteristics  and  to  correlate  with 
loss  in  stability  margin.  The  use  of  total  pressure  and 
total  temperature  has  been  shown  through  experience  and  is 
accepted  in  practice  to  be  the  most  practical  distortion 
pattern  description  method  (see  Bibliography) . The 
Society  of  Automotive  Engineers  Aerospace  Recommended 
Practice  on  Gas  Turbine  Engine  Inlet  Flow  Distortion 
Methodology  [12]  specifies  total  pressure  as  the  parameter 
for  distortion  pattern  description.  Total  pressure  and 
total  temperature  are  used  because  they  are  much  more 
readily  measured  than  velocity  distribution,  static  pressure 
and  flow  direction.  Also,  when  a distorted  flow  field 
approaches  an  engine,  the  interaction  of  the  engine  with  the 
distortion  pattern  produces  a specific  static  pressure. 
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velocity  and  flow  angularity  pattern  at  the  engine  entry 
(discussed  in  Chapter  IV) . Thus,  the  combination  of  a total 
pressure  or  total  temperature  distortion  pattern  with  an 
engine  at  a given  operating  condition  implicitly  specifies 
a velocity,  static  pressure  and  flow  direction  distribution 
at  the  compressor  entry  and  throughout  the  compressor. 
Measurement  of  velocity  and  flow  direction  at  the  compressor 
entry  is  sometimes  suggested  as  a more  physically  meaningful 
set  of  measurements  for  distortion  pattern  definition 
because  rotor  blade  angle  of  attack  is  related  to  the 
velocity  vectors  at  the  rotor  blade  leading  edge.  However, 
compressor  instability  caused  by  distortion  can  (and  often 
does)  occur  because  of  stall  in  any  stage  of  the  compressor. 
Therefore,  unless  detailed  velocity  and  flow  direction 
measurements  are  made  at  each  and  every  compressor  stage 
entry,  no  real  additional  information  is  provided  by  the 
more  difficult  velocity  and  flow  direction  measurements. 
Further,  the  total  pressure  and  total  temperature  distri- 
butions produced  by  an  aircraft  inlet,  if  measured  approxi- 
■Tiately  an  engine  diameter  ahead  of  the  compressor  entry 
plane,  will  bo  the  same  with  or  without  an  engine  installed 
behind  the  inlet.  Therefore,  distortion  pattern  information 
from  scale  model  and  full-scale  inlet  tests  may  be  obtained 
without  the  engine  installed.  This  could  not  be  done  if 
velocity  static  pressure  and  flow  direction  were  chosen  as 
the  distortion  pattern  definition  parameters.  For  these 
reasons,  total  pressure  and  total  temperature  variations 
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will  also  be  used  for  disturbance  description  for  the 
analytical  work  in  this  study. 

Most  of  the  influences  shown  in  Fig.  4 were 
determined  from  simulations  of  actual  flight  in  altitude 
test  cells  or  wind  tunnels  using  artificial  metnods  of  dis- 
turbance generation.  These  disturbances  are,  however, 
representative  of  actual  disturbances  experienced  in  flight. 
The  influence  on  stability  in  actual  flight  is,  however, 
more  difficult  to  measure;  thus,  the  ground  test  sources  are 
used.  More  detailed  information  on  the  generation  of  these 
disturbances  in  flight  is  given  in  [13  through  16] . 

III.  PREVIOUS  ANALYSES  OF  COMPRESSOR  STABILITY 
WITH  UNDISTURBED  FLOW 

Early  analysis  of  compressor  instability  dealt 
with  "surging"  described  as  "a  violent  fluctuation  in 
delivery"  by  Whittle  in  an  early  work  (1931)  on  the  turbo- 
compressor as  a supercharger  [17] . In  an  even  earlier  work 
by  Kearton  [18]  and  in  later  works  by  Den  Hartog  [19]  and 
Ilorlock  [20],  analysis  of  surge  was  performed  using  the 
shapes  of  the  pressure-ratio-mass  flow  curves  of  the  com- 
pressor and  the  throttle  (back  pressuring  device)  to  locate 
stable  and  unstable  combinations  of  the  compressor  and 
throttle  and  to  qualitatively  describe  surging.  (Figure  5 
is  a simple  illustration  of  this  analysis.)  With  the 
throttle  set  to  an  open  condition,  the  match  point.  A,  is  a 
stable  operating  condition.  This  may  be  deduced  by 
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imagining  that  the  flow  is  momentarily  perturbed  upward.  If 
this  should  happen,  the  upward  perturbation  will  cause  a 
reduction  in  the  pressure  exiting  from  the  fan,  and  an 
increase  in  the  resistance  to  flow  from  the  throttle.  Both 
of  these  effects  tend  to  cause  flow  to  decrease,  thus 
countering  the  original  perturbation  and  restoring  equi- 
librium. A similar  argument  applies  for  a downward 
perturbation  of  flow  rate.  If,  however,  the  match  point  is 
point  C,  instability  will  occur.  Once  more,  imagine  an 
upward  perturbation  of  flow.  This  will  result  in  an 
increase  in  compressor  discharge  pressure  which  is  greater 
than  the  increase  in  resistance  to  flow;  thus,  the  upward 
perturbation  in  flow  will  be  reinforced  causing  the  flow  to 
continue  to  increase.  Thus,  point  C is  unstable.  This 
analysis  may  be  generalized  to  state  that  the  compressor- 
throttle  combination  will  be  unstable  in  regions  where  the 
slope  of  the  compressor  characteristic  curve  is  greater  than 
the  slope  of  the  throttle  loss  characteristic  curve. 

Applying  this  rule,  if  operation  is  attempted  near  point  C, 
and  the  flow  is  perturbed  upward,  it  will  continue  to 
increase  until  it  reaches  point  B where  the  flow  increase  is 
countered  (that  is,  the  compressor  characteristic  slope 
becomes  less  than  the  throttle  curve  slope)  and,  ignoring 
dynamic  effects,  is  caused  to  begin  to  decrease.  Once  the 
flow  starts  down,  however,  it  will  continue  to  decrease 
until  it  reaches  a region  near  point  D where  the  slope  of 
the  compressor  characteristic  is  once  more  less  than  that  of 
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the  loss  characteristic.  Thus,  a limit  cycle,  qualitatively 
describing  surge  based  on  characteristic  curve  nonlinearity, 
but  ignoring  dynamic  effects,  is  constructed. 

This  simple  explanation  provided  a reasonable 
qualitative  explanation  of  the  surge  instability,  but  in 
actual  use  tended  to  predict  a value  of  flow  rate  at  which 
surge  would  be  encountered  which  was  often  lower  than  that 
actually  experienced  [21] . it  later  became  common  practice 
to  assume  that  instability  would  occur  at  the  peak  com- 
pressor ratio;  that  is,  at  the  first  point  where  the 
compressor  characteristic  slope  became  positive.  This 
assumption  was  reasonably  well  borne  out  for  high  corrected 
rotor  speeds  by  experiments  [21,  22]. 

While  the  theories  described  to  this  point  were 
based  principally  on  observations  made  on  centrifugal  com- 
pressors, they  are  in  principle  applicable  to  axial  flow 
units.  (For  example,  the  compressor  characteristic  curve  of 
^^9-  5 in  no  way  specifies  the  type  of  compressor  which 
produced  it.)  Pearson  and  Bowmer  (1949)  [21]  specifically 
addressed  the  problem  of  stability  in  axial  flow  compressors. 
Pearson  and  Bowmer  observed  the  discrepancy  which  existed 
between  compressor  experiments  and  the  simple  nonlinear 
characteristic  curve  theory  and  noted  that  the  nonlinear 
theory  ignored  the  dynamic  effects  caused  by  volumetric 
capacitance,  fluid  column  inertia,  and  flow  resistance  in 
the  compressor  and  ducting.  Pearson  and  Bowmer  approacned 
the  problem  by  deriving  an  approximate  set  of  dynamic 
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(time-dependent)  equations  for  a twelve-stage  axial  flow 
compressor  and  related  the  compressor  to  an  analogous 
electrical  circuit.  The  compressor  was  stacked  from  assumed 
stage  characteristics  to  form  a representative  multi-stage 
unit.  It  was  then  assumed  that  instability  would  occur  at 
any  condition  where  the  impedance  of  any  given  part  of  the 
circuit  representing  a stage  was  zero.  For  positive 
impedance,  disturbances  in  the  circuit  would  be  damped;  for 
negative  impedance  the  disturbances  would  be  amplified,  thus 
causing  instability. 

Howell  (1964)  [23]  extended  the  work  of  Pearson 

and  Bowmer  by  applying  the  same  technique  of  drawing  an 
electrical  analogy  to  the  compressor,  but  used  a more 
detailed  stability  analysis  on  the  resulting  electrical  net- 
work. Howell  also  applied  the  analysis  to  an  actual 
compressor  and  compared  those  results  to  the  experimentally 
determined  stability  limit.  Howell's  analysis  produced  a 
reasonably  good  prediction  of  the  stability  limit,  but  it 
was  somewhat  below  the  experimentally  observed  limit  at 
design  rotor  speed. 

Kimzey  and  Couch  (1969)  developed  a stability 
criterion  for  predicting  the  position  of  the  stability  limit 
directly  from  an  analysis  of  a simple  set  of  dynamic 
equations  representing  the  volumetric  capacitance  and 
inertial  characteristics  of  the  compressor  on  a stage-by- 
stage  manner  [24] . This  stability  criterion  was  based  on 
the  assumption  that  the  compressor  stages  work  in  pairs. 
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and  if  any  "pair"  becomes  an  unstable  combination,  then 
instability  of  the  entire  compressor  will  result.  (Any 
given  stage  in  the  compressor  is  a member  of  two  pairs;  one 
is  formed  in  conjunction  with  the  preceding  stage  and  one 
with  the  following  stage;  for  the  last  stage,  the  ducting 
and  throttle  form  one  pair;  for  the  first  stage,  the  up- 
stream ducting  forms  a pair  with  the  first  stage.)  The 
criterion  was  applied  to  an  actual  compressor  and  compared 
to  experimental  data.  Reasonably  good  agreement  was  shown. 

Daniele,  Blaha,  and  Seldner  (1974)  125]  applied 
more  sophisticated  stability  criteria  to  an  approximate  set 
of  linearized  dynamic  equations  representing  volumetric 
capacitance,  inertial,  and  energy  effects.  These  criteria 
were  applied  to  a compressor  for  which  detailed  experimental 
data  were  available  to  predict  the  steady-state  stability 
limit  line.  Very  good  agreement  with  the  experimental  data 
was  achieved. 

IV.  PREVIOUS  ANALYSES  OF  COMPRESSOR  STABILITY 
WITH  DISTURBED  FLOW 

The  work  described  above  was  principally  directed 
at  the  prediction  of  the  compressor  stability  limit  with 
steady  and  uniform  flow  entering  the  compressor.  Work  has 
also  been  performed  in  an  effort  to  predict  the  influence  of 
external  disturbances  on  the  stability  limit.  An  excellent 
and  comprehensive  review  of  the  work  accomplished  to  date  in 
this  area  is  given  in  [16] . 
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Steady  Nonuniform  (Di:;t.orto(i)  I’low 

Notable  early  work  dealing  with  the  influence  of 
distortion  at  the  compressor  entry  was  accomplished  by 
Alford  (1957)  [26]  . Alford  treated  the  problem  of  circumfer- 
ential total  pressure  distortion  with  a penetration  theory 
which  postulated  that  the  influence  of  circumferential 
distortion  on  stability  margin  was  proportional  to  the  depth 
of  penetration  of  the  distorted  pattern  into  the  compressor . 
The  distortion  pattern  attenuation  characteristics  were 
computed  in  an  approximate  manner,  thereby  deriving  the  form 

I 

of  an  index  which  correlated  distortion  pattern  character- 
istics to  loss  in  stability  margin.  Alford  concluded  that 
the  loss  in  margin  depended  on  the  magnitude  of  the  total 
pressure  deficit  in  the  low  pressure  region,  on  the  area  of 
the  compressor  entry  annulus  which  had  total  pressure  below 
the  average  value,  and  on  the  angular  extent  of  the  low 
pressure  region. 

Pearson  (1963)  [27]  described  the  parallel  com- 
pressor theory  of  circumferential  distortion.  The  theory 
divides  the  compressor  annulus  into  arc  sectors,  each  sector 
being  viewed  as  an  independent  compressor  running  in 
parallel  with  the  other  sectors.  All  sectors  are  assumed  to 
discharge  to  a common  static  pressure.  The  pressure  ratio- 
flow  rate  characteristic  curve  of  each  sector  was  assumed  to 
be  unchanged  by  the  distortion.  Crossflows  between  the 
adjacent  parallel  compressor  sectors  were  ignored.  Each 
sector  will  have  a different  inlet  total  pressure  because  of 
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the  circumferential  distortion  existing  at  the  compressor 
entry  plane.  In  its  simplest  form  of  application,  the 
entire  compressor  is  assumed  to  roach  its  stability  limit 
when  the  stability  limit  of  any  sector  is  encountered.  In 
practice,  correlating  indexes,  such  as  described  in  [28]  are 
developed  which  relate  the  actual  loss  in  stability  margin 
to  the  size  of  the  arc  sector  chosen  and  the  size  of  the 
deficiency  of  total  pressure  in  the  lowest  pressure  sector. 
Total  temperature  distortion  has  also  been  treated  using 
parallel  compressor  theory  [28]  . 

Goethert  and  Kimzey  [29]  presented  an  analog 
computer  model  for  calculating  the  effect  of  circumferential 
distortion  which  was  principally  based  on  the  parallel  com- 
pressor theory , but  did  make  provisions  for  approximating 
crossflow  effects  between  sectors.  Fett  (1969)  [24] 
combined  parallel  compressor  theory  with  the  Hurwitz 
stability  criteria  to  provide  a method  of  directly  pre- 

the  loss  of  stability  margin  caused  by  circumfer- 
ential distortion  for  simple  two-  and  three-stage 
compressors . 

Comparisons  between  the  degree  of  distortion 
predicted  to  cause  instability  by  the  parallel  compressor 
theory  in  its  simplest  form  and  that  observed  experimentally 
indicated  that  as  the  angular  extent  of  the  low  pressure 
region  was  reduced  from  180  degrees,  the  compressor  would 
actually  become  more  tolerant  to  the  distortion  than  pre- 
dicted. 
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Ucid  (1969)  [30]  , and  Calogeras,  Mchalic,  and 

Burstadt  (1971)  [31]  showed  experimentally  that  as  the 
angular  extent  of  the  low  pressure  region  was  increased  from 
0 to  60  degrees,  the  loss  in  stability  margin  increased 
drastically,  but  the  loss  was  nearly  constant  between  90  and 
180  degrees.  The  experiment  also  showed  that  somewhere 
between  60  and  90  degrees,  the  prediction  of  loss  of 
stability  margin  from  parallel  compressor  theory  was  more 
reasonable.  From  these  observations,  a critical  angle  was 
defined  as  the  minimum  distortion  angle  which  would  allow 
prediction  of  stability  margin  loss  by  simple  parallel  com- 
pressor theory. 

Goethert  and  Reddy  (1971)  [32]  arrived  at  a 

critical  angle  of  approximately  60  degrees  based  on  a 
theoretical  treatment  of  a cascade  of  blades  in  a phase- 
shifted  oscillatory  flow  field  such  as  would  be  produced  by 
a compressor  rotor  moving  through  a circumferentially  dis- 
torted flow  field.  In  the  unsteady  cascade  flow  treatment, 
it  was  shown  that  an  airfoil  (or  cascade  of  airfoils)  talces 
a certain  amount  of  time  to  change  its  lift  in  response  to 
rapid  changes  of  angle  of  attack.  It  was  also  shown  that 
the  change  in  the  lift  from  the  steady-state  value  during  a 
rapid  oscillation  of  the  approach  flow  direction  diminishes 
as  the  frequency  of  oscillation  is  increased.  Thus,  cir- 
cumferential distortions  with  low  pressure  regions  which 
cover  only  small  angular  extents  will  be  recognized  by  a 
rotor  blade  passing  through  the  region  as  a very  rapid 
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change,  and  the  change  in  lift  will  be  small.  As  the 
angular  extent  of  the  distortion  is  increased,  however,  the 
rotor  blade  will  have  an  increasingly  longer  time  to  adjust 
to  the  higher  angles  of  attack  existing  in  the  low  pressure 
region. 

The  portion  of  the  compressor  in  the  low  pressure 
region  will  therefore  operate  at  angles  of  attack  nearer  the 
stall  value,  thus  reducing  the  stability  margin  of  the 
compressor. 

The  analysis  of  Goethert  and  Reddy  was  particu- 
larly useful  in  that  the  difference  in  the  response  of  a 
cascade  compared  to  that  of  an  isolated  airfoil  was  clearly 
made.  Also,  the  difference  between  oscillatory  flow 
approaching  a stationary  cascade  and  stationary  flow 
approaching  an  oscillating  cascade  was  made.  The  influence 
of  unsteady  airfoil  aerodynamics  on  compressor  behavior  in  a 
distorted  flow  field  was  recognized  and  used  by  other 
investigators  also  [33,  34].  Extensions  have  also  been  made 
into  unsteady  airflow  aerodynamics  of  stalled  flow  around 
an  airfoil  by  Carta  (1973)  [35]  and  Melick  (1973)  [36] . 

Uniform  Time-Varying  Entry  Flow 

The  effect  of  uniform  time-dependent  disturbances 
on  the  stability  of  compressors  has  also  been  investigated. 
Gabriel,  Wallner,  and  Lubick  (1957)  [37]  studied  the  effects 
of  uniform  fluctuations  of  compressor  inlet  total  pressure 
on  compressor  stability.  The  influence  of  rapid  reunps  in 
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compressor  inlet  total  temperature  was  also  studied. 

Analysis  of  the  influence  was  performed  by  solving  a set  of 
lumped  volume  unsteady  mass  conservation  equations  along 
with  ec[uations  representing  the  pressure-ratio-mass  flow 
characteristic  curves  of  several  compressor  stage  groups. 

The  characteristic  curves  of  the  stage  groups  were  obtained 
from  steady-state  tests  of  the  compressor  and  were  assumed 
to  be  valid  dynamically  in  this  analysis.  This  is  the 
so-called  quasi-steady-state  stage  characteristic  assumption. 
The  behavior  of  a fifteen-stage  turbojet  engine  compressor 
was  analyzed  by  dividing  the  stages  into  four  stage  groups. 
The  resulting  equations  were  solved  using  an  analog  computer. 
It  was  assumed  in  the  analysis  that  if  any  stage  group 
characteristic  curve  reached  its  peak  value,  that  is,  the 
characteristic  curve  slope  changes  from  negative  to  posi- 
tive, that  compressor  instability  would  occur.  The  work 
indicated  that  for  uniform,  sinusoidal  oscillations  of 
compressor  inlet  total  pressure,  the  amplitude  required  to 
cause  instability  decreased  as  the  oscillation  frequency 
increased.  Also,  it  was  found  that  upward  ramps  in  com- 
pressor entry  total  temperature  would  cause  compressor 
instability,  and  that  the  effects  become  increasingly 
severe  as  the  ramp  rate  is  increased.  The  calculated 
results  were  compared  to  experimental  data  and  showed 
reasonably  good  agreement  with  the  experimentally  determined 
stability  limit.  This  particular  analysis  was  limited  to 
compressor  rotor  speeds  near  the  design  value. 
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Later  investigators  extended  Gabriel,  Wallner,  and 
Lubick's  approach  for  calculating  the  effects  on  stability 
of  uniform,  time-dependent  compressor  inlet  disturbances. 
Kimzey  (1966)  [38]  improved  the  analysis  by  using  each  com- 
pressor stage  separately  (instead  of  combining  into  stage 
groups) , and  improved  upon  the  form  of  the  analog  computer 
solution.  Further  refinements,  including  the  approximation 
of  unsteady  momentum  effects,  unsteady  energy  effects, 
improvement  of  the  criteria  for  instability,  and  imple- 
mentation of  the  solution  on  the  digital  computer,  have  also 
been  accomplished  [24,  39,  40]. 

Time-Varying  Distortion 

Although  a great  deal  of  excellent  experimental 
work  has  recently  been  accomplished  on  time-varying 
distortion  [13  through  16] , relatively  little  has  been 
accomplished  to  provide  a method  of  directly  calculating  its 
effects  on  stability.  A quasi-steady  approach,  developed  by 
Plourde  and  Brimlow  (1968)  [7] , is  most  commonly  used  at  the 
present  time.  This  particular  approach  assumes  that  the 
compressor  reacts  to  time-variant  distortion  in  the  same 
manner  as  it  reacts  to  steady-state  distortion  if  the  time- 
variant  distortion  pattern  "dwells"  at  the  compressor  entry 
for  a time  period  approximately  equal  to  the  time  required 
for  the  compressor  rotor  to  make  one  revolution.  A dis- 
tortion index,  which  correlates  the  loss  in  stability  margin 
to  the  spatial  distortion  pattern  characteristics  in  the 
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same  manner  as  is  used  with  steady-state  distortion,  is  then 
calculated  as  a function  of  time.  Time-dependent  fluctua- 
tions in  the  compressor  entry  flow  field,  which  are  at  a 
frequency  higher  than  that  commensurate  with  the  pattern 
"dwell  time,"  are  filtered  out  prior  to  making  the  dis- 
tortion index  computations.  The  loss  in  stability  margin  is 
then  correlated  with  the  maximum  distortion  index  value 
observed  over  a given  time  period.  Although  Plourde  and 
Brimlow's  treatment  does  not  directly  account  for  some 
dynamic  processes  known  to  be  present  in  the  compressor  with 
time-varying  distortion,  reasonably  good  correlation  is 
achieved  with  the  method,  making  it  a most  valuable  engi- 
neering tool  for  practical  applications. 

Melick  (1973)  [36]  applied  unsteady  airfoil  theory 
to  the  analysis  of  the  time-variant  distortion  problem. 
Additional  detail  on  time-varying  distortion  analysis  is 
provided  in  [41,  42]. 

V.  OBJECTIVES  AND  APPROACH  OF  CURRENT  WORK 

As  previously  stated,  the  overall  objective  of 
this  study  is  to  develop  an  improved  analytical  method  with 
general  applicability  for  predicting  compressor  stability 
margin  loss  caused  by  the  aforementioned  external  dis- 
turbances. The  development  of  the  method  is  intended  to 
provide  insight  into  some  of  the  physical  processes  which 
lead  to  compressor  aerodynamic  instability. 
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The  work  described  herein  specifically  extends  the 
present  state  of  knowledge  through  the  following  develop- 
ments . 

Improved  One -Dimens ional , Time -Dependent  Model 

An  improved  model , suitable  for  steady-state 
stability  limit  computation  and  for  analysis  of  planar, 
transient  and  dynamic  disturbances  is  developed.  The  model 
includes  mass,  momentum,  and  energy  effects.  The  full  non- 
linear form  of  the  equations,  not  used  in  earlier  works,  is 
solved.  Unsteady  cascade  influences,  also  not  included  in 
works,  based  on  the  Goethert-Reddy  analysis  [32]  , 
are  included  in  the  analysis  and  their  influences  evaluated. 
The  Goethert-Reddy  unsteady  cascade  analysis  is  compared  to 
experimental  unsteady  compressor  blade  aerodynamic  data. 
Detailed  comparisons  of  model  computed  steady-state 
stability  limits  are  presented  for  three  different  com- 
pressors providing  an  indication  of  the  generality  of  the 
analysis  method.  Four  unsteady  cases:  compressor 

instability  caused  by  planar  oscillatory  inflow,  dynamic 
unstalled  response  of  a compressor  to  oscillating  inflow, 
dynamic  unstalled  response  to  oscillating  back  pressure,  and 
response  of  a compressor  to  rapid  inlet  temperature  ramps, 
are  computed.  Where  possible,  the  computed  results  are 
compared  to  experimental  results. 
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Extension  to  Multi-Dimensional  Flow  Case 

The  concepts  of  the  one-dimensional,  time- 
dependent  model  are  extended  to  three  dimensions  to  allow 
treatment  of  distorted  flows  into  the  compressor.  The 
compressor  is  divided  both  radially  and  circumferentially 
into  control  volumes,  and  crossflows  among  the  volumes 
caused  by  distortion  and  their  effects  on  stability  are 
computed.  The  influence  of  the  radial  work  variation  is 
approximated  and  included  for  the  radial  distortion  evalu- 
ations. Unsteady  cascade  effects  are  included.  Thus,  a 
model  capable  of  computing  the  influences  of  combined  radial 
and  circumferential  distortion  is  produced.  The  model  is 
used  to  compute  the  loss  in  stability  margin  for  a com- 
pressor with  steady-state,  combined  radial  and  circumferen- 
tial distortion,  and  comparisons  to  experimental  results  are 
made.  A time -dependent  distortion  case  is  computed. 
Parametric  influence  of  pure  radial  and  pure  circumferential 
pressure  distortion  are  analyzed.  Pure  circumferential 
temperature  distortion  is  also  analyzed. 

Comparisons  to  Experimental  Results 

Emphasis  was  placed  on  comparisons  of  model 
computations  to  experimental  results  to  maximize  under- 
standing of  the  physical  processes  and  to  test  the  validity 
of  the  computations.  Experimental  results  on  no  single 
compressor  were  sufficient  for  the  comparisons  desired; 
thus,  three  different  compressors  were  modeled.  They  are: 
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a four-stage  Allison  XC-1  compressor  [43] , an  eight-stage 
NACA  research  compressor,  the  NACA-8  [44,  45,  46],  and  the 
General  Electric  eight-stage  J85-13  engine  compressor  [47], 
Information  on  these  compressors  used  in  this  study  are 
given  in  Appendix  B.  The  NACA-8  and  General  Electric  J85-13 
compressors  were  modeled  one-dimensionally  and  the  models 
included  a portion  of  inlet  ducting  and  combustor  ducting 
for  more  realistic  time-dependent  boundary  conditions  for 
the  compressor.  The  Allison  XC-1  compressor  was  used  for 
one-dimensional,  steady-state  stall  analysis  and  for  the 
distortion  studies.  No  ducting  was  included  in  the  XC-1 
models. 


VI . ORGANIZATION 

The  development  of  the  one-dimensional,  time- 
dependent  model  is  presented  in  Chapter  II.  Chapter  III 
presents  application  of  the  one— dimensional  analysis  to 
example  problems . The  extension  of  the  model  to  three 
dimensions  is  given  in  Chapter  IV.  Chapter  V describes 
application  to  distortion  cases.  The  work  is  summarized 
and  recommendations  for  further  work  is  given  in  Chapter  VI. 
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CHAPTER  II 

DEVELOPMENT  OF  ONE-DIMENSIONAL, 

TIME-DEPENDENT  MODEL 

The  axial-flow  compressor  is  made  up  of  successive 
rows  of  airfoils,  a moving  row,  the  rotor,  followed  by  a 
stationary  row,  the  stator,  as  shown  in  Fig,  6.  A rotor 
followed  by  a stator  is  defined  as  a stage  for  this  study. 

A representative  absolute  velocity  and  static  pressure 
distribution  through  the  compressor  is  also  shown  in  Fig.  6. 
The  compressor  increases  the  static  pressure  and  density  of 
the  fluid  passing  through  it  by  the  dynamic  action  of  the 
rotor  which  imparts  kinetic  energy  to  the  fluid.  Depending 
on  the  specific  design,  the  static  pressure  through  the 
rotor  changes  by  varying  amounts.  The  function  of  the 
stator  row  is  to  diffuse  and  redirect  the  flow  to  convert  a 
portion  of  the  imparted  kinetic  energy  to  internal  energy. 
The  flow  is  further  diffused  as  it  passes  from  the  com- 
pressor through  the  compressor  diffuser  to  the  combustor. 
After  leaving  the  combustor  it  passes  through  a turbine 
nozzle  (if  on  an  engine)  or  a throttling  valve  (if  on  a 
compressor  test  rig)  which  is  usually  choked. 

In  the  formulation  of  the  one-dimensional  model,  a 
control  volume,  enclosing  the  fluid  in  the  compressor  is 
drawn,  as  shown  by  the  dashed  lines  in  Fig.  7.  For  the 
one-dimensional  analysis,  all  flow  properties  are  "bulked" 
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in  the  radial  and  circumferential  directions  and  variations 
with  the  axial  coordinate  and  with  time  only  are  considered. 

The  development  of  the  one-dimensional  model  pro- 
ceeds as  follows.  The  governing  equations  are  written.  The 
method  of  specifying  the  force  and  shaft  work  acting  on  the 
fluid  is  developed.  The  finite  difference  forms  of  the 
equations  used  for  numerical  solution  on  a digital  computer 
are  then  developed  and  the  method  of  solution  given. 

I . GOVERNING  EQUATIONS 

Figure  7a  illustrates  the  ducting-compressor- 
burner-turbine  vane  system  which  was  modeled.  Figure  7b 
shows  the  overall  control  volume  representing  the  system. 

The  forces  of  the  blading  and  walls  of  the  system  acting  on 
the  fluid  are  represented  by  an  effective  axial  time- 
dependent  force  distribution,  F(z,t).  Similarly,  shaft 
work  done  on  the  fluid  and  heat  added  to  the  fluid  are 
represented  by  the  time-dependent  distributions,  WS{z,t)  and 
Q(z,t) . Mass  removal  or  addition  (interstage  bleed,  for 
example)  are  represented  by  another  time-dependent  bleed 
flow  distribution  function,  WB(z,t).  All  bleed  flows  are 
assumed  to  cross  the  control  volume  boundary  normal  to  the 
axial  direction.  Time-dependent  boundary  conditions  are 
provided  by  specification  of  the  total  pressure  and  total 
temperature,  P(t)  and  T(t),  respectively,  at  the  system 
forward  boundary  and  by  specification  of  discharge  static 
pressure,  ^gxit^^^'  boundary.  For  choked  turbine 
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vanes,  a sufficiently  low  discharge  static  pressure  was 
specified  to  assure  choking.  In  that  case,  the  aft  boundary 
condition  is  unity  Mach  number.  If  unchoked,  specification 
of  static  pressure  is  the  aft  boundary  condition. 

The  governing  equations  are  derived  by  appli- 
cation of  time-dependent  mass,  momentum,  and  energy 
conservation  principles  to  the  elemental  volume  of  Fig.  7c. 


Mass 


Application  of  the  mass  conservation  principle 
(mass  is  neither  created  nor  destroyed)  to  the  elemental 
control  volume  yields 


W + dz  + WB  dz  + 

dZ 

mass  leaving 
control  volume 
per  unit  time 


3 (pA  dz) 

at 


time  rate  of 
increase  of  mass 
"stored”  in  the 
control  volume 


W . (1) 

- 

mass  entering 
control  volume 
per  unit  time 


Equation  (1)  may  be  reduced  to 


3 (pA)  _ ^ 

at  3z 


WB 


(2) 


Momentum 

The  momentiim  principle , the  summation  of  all 
forces  in  a given  direction,  is  equal  to  the  time  rate  of 
change  of  momentum  in  that  direction,  gives 
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F dz  + PS  A - 

PS  A + — ^ dz 

+ PS 

-A 

L J 

- 

- 

axial  forces  acting  on  control  volume 


momentum 

leaving 

control  volume 
per  unit  time 


WU 


momentum 
entering 
control  volume 
per  unit  time 


time  rate 
of  increase 
of  momentum 
"stored"  in  the 
control  volume 


(3) 


total  time  rate  of  change  of  momentum 


Equation  (3)  may  be  reduced  to  give 


3W 


(4) 


where,  the  impulse,  IMP,  is  defined  as 


IMP  = WU  + PS  A 


(5) 


Energy 


Energy  conservation  gives. 


enthalpy  leaving 
control  volume 
per  unit  time 


time  rate  of  increase 
of  energy  "stored"  in 
control  volume 


+ WS  dz  + Q dz 


H 

enthalpy 
entering 
control  volume 
per  unit  time 


> ^ 

shaft  work 
done  on 
fluid  in 
control  volume 


heat  added 
to  fluid  in 
control  volume 


(6) 
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Equation  (6)  may  be  reduced  to 


where 


and 


3(XA)  _ 

at 


-|2  + ws 

9z 


+ Q 


H = c WT  . 
P 


(7) 


(8) 


(9) 


State  and  Additional  Equations 

Additional  equations  required  include  the  perfect 
gas  equation  of  state. 


PS  = D R TS  , (10) 

calorically  perfect  internal  energy  and  enthalpy  relation- 
ships , 


e = c^  TS  + constant  , 

(11) 

h = Cp  TS  + constant  , 

(12) 

stagnation  (total)  state  temperature. 

2 

T = Cp  TS  + ^ , 

(13) 

and  stagnation  (total)  state  pressure. 

r T 

(14) 
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where 


c 


Also  required  is  the  Mach  number, 


(15) 


M--  , 


(16) 


where  the  acoustic  velocity,  a,  is 


a = /Y  R TS  = -)Jy  ^ , 


(17) 


and  the  Mach  number-static  pressure  form  of  the  energy  in 
the  control  volume. 


which  results  from  combination  of  Eqs.  (8),  (10),  (11), 
(15) , (16) , and  (17) . 

The  Mach  number-static  pressure  form  of  the 
Impulse,  IMP,  is  also  used, 

IMP  = PS  A [ 1 + yM^]  . 


(18) 


(19) 


This  equation  results  from  combination  of 


W = pAU  , 


(20) 


with  Eqs.  (5)  an6  (17). 
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II.  STACL  I-'ORCES  AND  SHAFT  WORK 


Steady-State  Stage  Characteris-tics 

Consider  the  element  of  fluid  contained  between 

two  blades  in  a rotor  as  shown  in  Fig.  8a.  Entering  the 

rotor,  the  tangential  velocity  is  Vp  , and  at  the  exit  it  is 

“l 

Vp  . Applying  the  conservation  of  angular  momentum 
principle  to  the  fluid  element  in  steady  flow  gives  the 
torque  on  the  element, 


i = W 


^1 


(21) 


Multiplying  Eq.  (21)  by  rotor  speed,  u,  yields  the  time  rate 
of  doing  work  on  the  fluid,  which  is  equal  to  the  total 
enthalpy  change  across  the  rotor.  Thus, 

''Oj]  • <”) 


Replacing  the  difference  between  radius -velocity  products 
using  an  effective  radius  yields. 


where 


(23) 


r 


V. 


(24) 
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Wrl'ting  wheel  speed  as 


Uwh  = “^  (25) 

and 

^02  = ''''6  «26) 

gives 

T_  - T,  U , 

'^''e  • <2’> 

Applying  the  linear  momentuin  equation  for  steady 
flow  to  the  element  of  Fig.  fla  gives 

Fg  =W  AVg  = U AVg  , (28) 

where  Fg  is  the  force  being  applied  to  the  fluid  by  the 
blading  in  the  compressor. 

The  force,  Pg , can  also  be  expressed  in  terms  of 
blade  lift  and  drag,  as  indicated  in  Fig.  8b, 

Fg  = L cos(a+A)  + D sin(o+A)  , (29) 

where  a is  the  blade  angle  of  attac)<,  determined  by  the 
velocity  triangle  of  Fig.  8b,  and  X is  the  blade  chord-to- 
axial  direction  angle,  or,  the  cascade  stagger  angle. 
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The  lift  and  drag  may  be  expressed  as  coef- 
ficients. 


and 


Pi  V , A - 
1 rel  ref 


(30) 


^dT 


(31) 


t Pi  '^rel  *ref 


Combining  Eqs.  (27)  through  (31)  produces 


T - T 
2 


_ 1 *ref  Hiiii  rel  T 

2Cp  Aj  J H [ 


oos(a+A)  + sin(a+A) 


(32) 


Returning  once  more  to  Fig.  3b,  define  the  stage  flow 
coefficient  as 

^ . (33) 

''e 

From  Fig.  8b,  the  flow  coefficient,  4>i  is  related 
to  the  blade  angle  of  attack,  a,  by  the  stagger  angle,  A, 
which  is  a constant  for  any  given  blade  row, 

= cot  (A+a)  . (34) 

Thus,  $ serves  the  same  purpose  in  the  compressor  as  a does 
for  aircraft  wing  aerodynamics.  It  decreases  with 
increasing  a;  therefore,  small  <>  means  large  a,  and  vice 
versa.  Also  from  Fig.  8b, 
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u 


rel 


= cos  ( A+fx) 


(35) 


Combining  £qs.  (32)  through  (35)  yields 


1 A , 
1 ref 


wh 


(36) 


The  left-hand  side  of  Eq.  (36)  is  defined  as  the  stage 
loading  parameter,  (following  Horloc)t  [20]),  or  simply  as 
the  stage  temperature  coefficient  (more  common  in  American 
worlcs  [40]).  Thus, 


c 

JS. 


(TR-1) 


(37) 


where  TR  is  the  stage  total  temperature  ratio. 

If  the  swirl  component  of  tangential  velocity  is 

small. 


and 


d>  = 


U. 


U 

'wh 


(38) 


(39) 
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At  given  Mach  and  Reynolds  numbers,  and  are  functions 
of  a;  thus,  ((>,  only.  Therefore, 


T 


Re,M 


(40) 


In  principle,  theoretical  or  experimental  lift  and 
drag  coefficients  could  be  used,  along  with  stage  geometry, 
(A,  to  compute  the  stage  temperature  coefficients 

following  Eqs.  (36)  or  (39).  Corrections  for  blade  row 
interference  causing  losses  and  secondary  flows  would  be 
required.  Another  method  is  to  measure  the  stage  total 
temperature,  flow  rate,  total  pressure,  and  flow  angularity 
at  the  stage  entry  and  exit  and  directly  compute  euid  $ 
from  Eqs.  (37)  and  (38).  The  stage  characteristics  used  in 
this  study  were  obtained  in  that  manner  and  are  summarized 
in  Appendix  B. 

The  temperature  ratio  of  Eq.  (37)  for  an  ideal 
compressor  stage  is  the  isentropic  temperature  ratio  defined 
from  stage  entry  and  exit  stagnation  pressure.  Using  this 
ratio,  another  stage  parameter,  the  stage  pressure  coef- 

p 

ficient,  i|j  , is  produced. 


P 


c 

P 


] 


(41) 
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The  pressure  coefficient  is  also  a function  of  the  flow 
coefficient,  ((i,  in  a manner  similar  to  the  temperature 
coefficient.  The  stage  adiabatic  efficiency  is 


'ad 


Ah. 

isen tropic 
^ actual 


I±1 

PR  ^ - 1 
TR=I 


Thus , 


A set  of  curves, 


1 = ^ 

,|,T 


m fTi 

((P) 


P P 

r =r  (<J>) 


> . 


'^ad  = ^ad 


(42) 


(43) 


(44) 


is  referred  to  as  stage  characteristics,  and  along  with  flow 
direction. 


6 = tan 


(45) 


or  an  equivalent  set  of  flow  direction  information,  fully 
defines  a stage's  performance  from  a one-dimensional  stand- 
point. Note  that  the  set  of  Eq.  (44)  is  redundant;  any  two 
of  the  three  variables  are  sufficient  to  compute  the  third. 

Various  forms  of  the  stage  characteristics  are 
used.  An  alternate  form  used  in  NACA  compressor  research. 
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and  used  for  the  stage  characteristics  of  the  NACA-8  com- 
pressor treated  in  this  study,  is  the  equivalent  temperature 
and  pressure  ratio  form.  They  are  defined  as 


and  is  defined  as  in  Eq.  (43) . 
Because 


nkl 

*^p^  design 


(46) 


(47) 


is  a constant  for  a given  compressor,  the  set  remains  a 
function  of  as  in  Eg.  (44),  and  contains  the  same 
information. 

Figure  9 is  an  example  set  of  stage  character- 
istics for  the  first  stage  of  the  NACA-8  compressor  [46] . 

To  the  right  of  the  equivalent  pressure  ratio  peak,  the 
stage  is  unstalled  (i.e.,  high  ^ and  low  a).  To  the  left, 
the  stage  is  stalled.  Rapid  decrease  of  stage  efficiency  at 
low  (])  (high  a)  where  the  blades  are  stalled,  is  evident. 

The  data  of  Fig.  9 is  seen  to  generalize  well  with  the  flow 
coefficient,  computed  as  in  Eg.  (38)  with  the  swirl 
component  of  velocity  ignored.  This  implies  that  for  this 
set  of  data,  the  swirl  velocity  component  was  either  small 
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relative  to  wheel  speed,  or,  changed  little  as  a fraction  of 
wheel  speed  over  the  range  of  compressor  operation  con- 
sidered. Steady-state  stage  characteristics  for  the  three 
compressors  treated  in  this  study  were  obtained  from  [43, 

46,  47],  and  are  shown  in  Appendix  B. 

The  foregoing  discussion  pertains  primarily  to  the 
compressor  rotor  which  imparts  the  shaft  work  to  the  fluid. 
No  stagnation  temperature  change  occurs  across  the  stators 
(with  adiabatic  boundaries)  and  the  stagnation  pressure 
losses  of  the  stator  may  be  combined  with  the  change  across 
the  rotor  to  give  a net  change.  Therefore,  the  character- 
istics of  Eq.  (44)  may  be  attributed  to  the  overall  stage. 

Unsteady  Cascade  Effects 

When  the  angle  of  attack  of  the  flow  into  the 
compressor  stage  varies  slowly,  the  lift  coefficient  follows 
the  steady-state  Cj^-a  curve  very  closely,  as  illustrated  in 
Fi?*  10 ■ Similarly,  the  steady-state  ifi— characteristics 
are  followed.  When  a rapidly  changes,  however,  the  lift 
cannot  instantaneously  follow  because  a certain  amount  of 
time  is  required  for  the  flow  around  the  blade  to  adjust  to 
the  new  condition.  Lag  loops  or  hysteresis  in  the  C.-a  and 
curves  appear  (Fig.  10)  . 

Goethert  and  Reddy  analyzed  the  unsteady  cascade 
effects  [32]  . They  analyzed  unsteady  incompressible 
potential  flow  through  a stationary  cascade  of  thin,  flat 
plate  foils  with  small  oscillations  in  the  axial  velocity 
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producing  angle-of-attack  oscillations.  The  time-averaged 
angle  of  attack  was  zero.  Their  results  are  summarized  in 
Fig.  11.  The  cascade  stagger  angle  is  A;  the  blade  chord 
length,  C;  and  spacing  is  a.  Rotor  speed,  u,  was  constant 
and  was  zero.  The  axial  velocity,  U,  was  oscillated  at 
frequency  f,  causing  an  inverse  oscillation  of  a.  This 
produced  an  oscillation  in  the  lift  coefficient.  The 
results  were  presented  in  terms  of  maximum  amplitude  ratio 
at  a given  frequency,  divided  by  the  lift  coefficient 
amplitude  at  f = 0 (quasi-steady-state) . The  resulting  lift 
ratio  is  presented  as  a function  of  reduced  frequency,  k, 
for  four  values  of  cascade  spacing,  a/C.  The  case  of 
a/C  = “,  a single,  isolated  airfoil,  is  shown  for  reference. 
Also  treated  for  reference  was  the  case  of  oscillating  foils 
in  a steady  stream.  The  Goethert-Reddy  work  showed  signifi- 
cant differences  to  exist  between  the  oscillating  foils, 
stationary  stream,  and  the  converse  case.  Their  work  also 
showed  the  unsteady  effects  for  a cascade  of  finite  spacing 
to  be  much  less  severe  than  that  of  an  isolated  airfoil. 

The  unsteady  influence  was  shown  to  increase  with  increasing 
spacing  (a/C)  and  stagger  angle,  1. 

The  lift  coefficient  ratio  variation  with  reduced 
frequency  is  intuitively  correct  at  k = 0 and  as  k 
Following  Goethert  and  Reddy  [32]  the  reduced  frequency,  k, 
may  be  viewed  as  an  indication  of  the  ratio  of  flow  passage 
time  to  the  disturbance  stay  time  of  the  forced  oscillation. 


37 


AEOC-TR-77-SO 


The  passage  time  is  defined  as  the  time  required  for  an  air 
psj^ticle  to  pass  through  the  blade  row  channel.  Or, 


At 


passage 


(48) 


The  disturbance  stay  time  is  defined  as  the  time  that  any 
one  blade  is  exposed  to  the  positive  disturbance  velocity  of 
one  velocity  oscillation,  that  is. 


Therefore, 


At 


^ 1 

disturbance  7F  ’ 


(49) 


At 


^^passage 

disturbance 


(50) 


Alternately,  k may  be  viewed  as  being  proportional 
to  the  ratio  of  the  blade  chord  length  to  the  disturbance 
wave  length. 


chord  length  _ C _ Cf  k . 

wave  length  '^rei  ~ '^rel 

If  k is  near  zero  (i.e.,  the  disturbance  time  period  is 
large  compared  to  the  passage  time) , the  blade  lift  will 
easily  follow  the  disturbance.  Similarly,  the  chord  length 
is  small  relative  to  the  disturbance  wave  length;  thus,  the 
lift  can  easily  respond.  If  k -^  <»,  the  disturbance  time  is 
small  compared  to  the  passage  time  and  many  disturbance  wave 
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lengths  reside  within  one  chord  length  at  any  instant  in 
time.  The  disturbance  is  then  "averaged-out"  by  the  blade 
yielding  an  instantaneous  lift  value  of  zero  due  to  the 
disturbance. 

Bruce  and  Henderson  [48]  conducted  an  experi- 
mental investigation  of  the  unsteady  force  coefficients 
produced  by  a compressor  rotor  in  unsteady  flow.  A single 
stage,  low-speed  compressor  was  fitted  with  a rotor  blade 
having  a segment  of  the  blade  on  a small  strain— gaged  force 
balance.  This  installation  allowed  a direct  measurement  of 
the  unsteady  forces  acting  on  that  blade  segment  as  the 
rotor  passed  through  a distorted  flow  field  produced  by  a 
distortion  screen.  A portion  of  their  measurements  was 
taken  with  stagger  angle  and  cascade  spacing  near  the  values 
used  by  Goethert  and  Reddy.  Figure  12  shows  a comparison  of 
the  Goethert-Reddy  analysis  with  the  Bruce-Henderson  experi- 
mental results.  The  trends  with  reduced  frequency  are  in 
good  agreement.  The  experimental  results  generally  were 
from  10  to  20%  higher  than  the  analytical  results.  Con- 
sidering the  difficulties  and  assumptions  involved  in  both 
the  analytical  and  experimental  evaluations,  this  level  of 
agreement  is  considered  reasonable  at  this  time.  Agreement 
among  other  theoretical  treatments  and  experimental  results 
surveyed  was  usually  of  this  order  or  worse.  (See  compari- 
sons made  in  [48]  for  example.) 

The  analysis  of  Goethert  and  Reddy  was  performed 
for  an  unloaded  cascade,  that  is,  the  mean  angle  of  attack 
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was  zero.  Figure  13  shows  the  variation  of  lift  coefficient 
ratio  with  mean  angle  of  attack  from  the  unsteady  cascade 
work  of  Ostdiek  [49],  and  Bruce  and  Henderson's  single-stage 
compressor  [48J  . The  dynamic  lift  coefficient  ratio  magni- 
tudes were  normalized  by  the  value  of  a mean  angle  of  attack 
of  eight  degrees.  The  measurements  were  made  over  a range 
of  reduced  frequency  of  0.002  to  2.9,  and  from  0 to  12 
degrees  mean  angle  of  attack.  The  results  are  all  within 
±10%  and  no  generalized  trend  could  be  identified. 

Unsteady  Stage  Characteristic  Correction 

The  Goethert-Reddy  analysis  may  be  used  to  con- 
struct a correction  to  the  steady-state  stage  character- 
istics. To  facilitate  application  of  the  correction,  the 
dynamic  lift  coefficient  ratio  was  represented  by  a simple, 
f j-J^st-order  differential  equation,  as  indicated  in  Fig.  14. 
The  equation. 


(52) 


represents  the  time-dependent  lift  coefficient,  C^,  in 
response  to  the  quasi-steady-state  lift  coefficient,  C.  , 

S S 

which  would  be  produced  if  the  lift  exactly  followed  the 
steady-state  Cj.-a  relationship. 

The  blade  time  constant,  t,  depends  on  the  chord 
length  and  reduced  frequency. 
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T 


£ C 


V 


rel 


(53) 


The  constant,  e,  depends  on  cascade  geometry  (i.e.,  a/C  and 
A).  For  a/C  = 1.0,  and  A = 50  degrees,  e = 0.3106.  The 
value  of  E is  found  by  fitting  the  approximating  differ- 
ential equation,  Eq.  (52),  to  the  response  indicated  by  the 
Goethert-Reddy  analysis . 

For  sinusoidal  variation  of  angle  of  attaclc,  the 
amplitude  response  of  Eq.  (52)  is 


max 


'£  ss 


/ — A . 

y 1 + (E 


(54) 


The  approximation  of  Eq.  (52)  is  compared  to  the  Goethert- 
Reddy  solution  in  Fig.  14,  and  fidelity  -is  excellent  up  to 
k = 8. 

The  cascade  airfoil  lift  coefficient  is  related  to 

the  stage  characteristics  through  Eq.  (39)  . Consider 

oscillations  of  the  flow  into  a stage  around  a given  stage 

characteristic  operating  point  0^)  • Let  0 represent 

T P 

exther  0 or  ij;  . Expanding  Eq.  (39)  in  a Taylor  series, 
neglecting  the  drag  coefficient  term  and  retaining  only  the 
first-order  terms  of  the  series  gives 


0-0, 


1 \ef 


2 A, 


(55) 
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Near  the  point,  (lU^,  (p^) , 


Then,  fromEqs.  (55)  and  (56), 


(56) 


C 


Z 


(57) 


Substituting  Eq.  (57)  into  Eq.  (52)  produces. 


^ i 

dt  T T * 


(58) 


Thus,  the  stage  characteristics  may  be  corrected  for 
unsteady  cascade  airfoil  influences  using  the  same  first- 
order  differential  equation  algorithm  as  that  constructed 
for  the  unsteady  lift  coefficient. 

The  linearization  of  the  unsteady  correction 
implicit  in  Eqs.  (55)  and  (56)  is  justified  by:  1)  the 

Goethert-Reddy  analysis  is  itself  a linearized  analysis; 

2)  the  correction  magnitude  is  small  in  the  range  of  reduced 
frequency  normally  encountered,  as  will  be  demonstrated 
later.  The  correction  is  actually  proper  only  up  to  the 
point  (i.e.,  up  to  the  point  where  unsteady  viscous 
effects  become  predominant) . 
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Force  and  Shaft  Work  Distributions 

The  momentum  equation,  Eq.  (4),  requires  an  axial 
force  distribution  representing  the  force  of  the  compressor 
blading  and  casings  on  the  fluid.  The  function  is  con- 
structed from  the  experimentally  determined  steady-state 
stage  characteristics  corrected  for  unsteady  effects.  (Or, 
alternately,  the  force  may  be  computed  in  a quasi-steady- 
state  manner  from  the  steady-state  stage  characteristics  and 
then  directly  corrected  for  unsteady  effects  using  the 
dynamic  lift  coefficient  ratio  of  Fig.  14.)  The  shaft  work 
distribution  required  by  the  energy  equation,  Eq.  (7),  is 
similarly  constructed. 

The  exact  method  of  extracting  the  force  and  shaft 
work  distributions  depends  partly  on  the  finite  difference 
method  used  to  solve  the  governing  equations,  and  is 
described  in  the  next  section. 

Ill,  METHOD  OF  SOLUTION 
The  basic  governing  equations  are: 


5 (pA) 
3t 


(2) 


and 


|W  . . DUMPI.  ^ p , PS  |A 

3t  dZ  o2  ' 


(4) 


M . MS 

or  d z 


+ Q 


(7) 
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The  dependent  variables  to  be  integrated  for  are  pA,  W,  and 
XA.  The  area.  A,  was  preserved  in  the  derivative  to  allow 
for  cases  where  exit  area  is  rapidly  changed. 

Finite  Difference  Approximation,  One-Sided 
Difference  Form 

The  spatial  derivatives  were  approximated  in  two 
different  finite  difference  forms.  When  only  the  compressor 
was  considered  (no  ducting) , a simple,  one-sided  difference 
approximation  was  used.  The  control  volume  was  divided  into 
unequal  lengths.  Each  length  included  a stage  of  the  com- 
pressor as  shown  in  Fig.  15.  Figure  15  also  shows  the 
computation  plane  designations.  The  spatial  derivatives  at 
location  i were  approximated  by 


3W. 

1 _ 

~5F  ■■ 


w. 

1 


^i 


(59) 


3 (IMP)^ 
3z 


IMP. 


ip 


^ip 


IMP^ 


z . 

1 


(60) 


3A. 

- constant  , (61) 

(the  product  of  static  pressure  with  the  area-derivative  was 

combined  with  the  force  term,  F) , 

and 


3H.  H.  - H. 

1 _ 1 im 

oZ  Z . - 2 . 

1 im 


(62) 


44 


AEDC-TR-77-80 


Substitution  into  Eqs.  (2),  (4),  and  (7)  produces 
the  approximate  equations. 


3 {pA) . W.  - w. 

ip  _ 1 ip 


3t 


z . - z . 

ip  1 


(63) 


3W.  IMP.  - IMP. 

= h ^ 4-  p 

3t  z . - z . ^ 1 ' 

ip  1 


and 


(64) 


3 (XA) . H.  - H. 

^ ^ + WS.  + Q. 

dt  Z.  - Z.  1 

xp  1 


(65) 


Writing  Eqs.  (63),  (64),  and  (65)  from  i = 1 to 
i = in-1;  that  is,  from  the  forward  to  the  aft  control 
volume  boundaries,  results  in  3x(in-l)  first-order  differ- 
ential equations  to  be  solved  simultaneously. 

Force  and  Shaft  Work  Computation 

The  division  of  the  control  volume  into  stage 
lengths  was  made  because  the  forces  and  shaft  work  applied 
between  stage  entry  and  exit  are  obtainable  from  the  stage 
characteristics.  To  determine  the  force  and  shaft  work 
applied  to  each  element,  the  flow  coefficient  into  the 
element  is  calculated, 


(66) 
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The  Stage  characteristic,  represented  as  a polynomial  curve 
fit,  is  then  entered  and  and  (or  are  found. 
Corrections  for  unsteady  effects  are  made.  The  stage  stag- 
nation pressure  ratio  and  temperature  ratios  which  go  into 
the  force  and  shaft  work  computations  are  calculated  from 


and 


Shaft  work  is  computed  as 


ws^  = W^T^  J TR^  - 1 j . 

Similarly,  the  axial  force  is  given  by 

= IMP^  ^ ^ j ‘ 


(67) 


(68) 


(69) 


(70) 


IMPR^  is  the  impulse  function  ratio  across  the  stage, 
determined  from  the  stagnation  temperature  and  pressure 
ratios,  area  ratio,  and  flow  directions. 


(71) 
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The  variable,  ME,  is  the  equivalent  Mach  number  downstream 
of  the  stage  which  would  be  present  if  steady  flow  existed. 
(Recall  that  the  stage  characteristics  from  which  the  forces 
and  shaft  work  are  determined,  are  based  on  experimental 
steady-state  measurements) . The  value  of  ME  is  computed 
from  a polynomial  expression  approximating  the  inverse  of 
£q.  (72) , for  constant  R and  y , 


MFFE^ 


+ 


Y-1 

2 


(72) 


MFFE  is  the  downstream  equivalent  mass  flow  function. 


MFFE. 

1 


w.  Vt.  tr. 

1*3.  1 

PR.  A.  COS 
1 ip 


(73) 


defined  in  the  same  spirit  as  was  HE. 

The  flow  direction  angles,  0^  and  taken  with 

respect  to  the  axial  direction,  are  input  as  functions  of 
compressor  corrected  rotor  speed,  as  constants,  or,  are 
assumed  to  be  zero  if  no  flow  direction  information  is 
available . 


Solution  Procedure 

Figure  16  outlines  the  overall  digital  computer 
solution  procedure  followed.  The  problem  is  solved  as  an 
initial  condition  problem.  Following  the  flow  path  of 
Fig.  16,  at  time  equal  zero,  initial  values  of  the  dependent 
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variables  are  specified.  Entry  plane,  interior  plane  and 
exit  plane  conditions  along  with  forces  and  shaft  work  are 
computed  for  time  equal  zero  and  are  used  to  compute  the 
time  derivatives  of  the  dependent  variables . The  time 
derivatives  are  then  integrated  using  a fourth— order  Runge- 
Kutta  numerical  integration  method  [50]  for  the  dependent 
variable  values  at  the  next  step  in  time.  The  new  values 
are  then  used  to  compute  conditions  for  that  time  step  and 
the  solution  advances  in  time  until  the  desired  solution  is 
obtained . 

Some  additional  details  of  the  computation  pro- 
cedure are  as  follows. 

Entry  plane  conditions . At  the  forward  boundary, 
i = if  flow,  will  be  available  from  integration  of 
Eq.  <64).  Density,  and  the  energy  variable,  are 
not  available  and  must  be  supplied  by  the  boundary  condition. 
Stagnation  pressure,  , and  stagnation  temperature,  are 
supplied  as  boundary  conditions.  The  mass  flow  function, 


MFF^ 


w.  VI7 
-rpq- 


(74) 


is  computed.  Mach  number  is  related  to  the  mass  flow 
function  by 


(75) 
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The  inverse  form  of  Eq.  (75)  is  approximated  by  a polynomial 
curve-fit  of  as  a function  of  MFF^  for  constant  R and  y. 
Static  temperature,  static  pressure,  density,  and  the  energy 
variable  may  be  computed  knowing  Mach  number: 


(76) 


(77) 


and 


(78) 


(79) 


Interior  plane  values.  On  the  interior  planes, 
values  of  the  three  dependent  variables  are  known  from 
integration.  To  provide  the  other  variables  necessary  for 
the  solution,  the  following  calculations  are  made  at  each 
plane : 


(80) 


(81) 


PS^  = R TS^  , 


(82) 
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and 


T. 

1 


TS^ 


r 


Y 


^Jy  ^ Ts.  , 


(83) 


(84) 

(85) 


U. 


IMP.  = 


1 - PSi  Ai[l  + T M?]  , 


(86) 

(87) 


H. 


0. 


c W.  T. 
p 1 1 


(88) 


Exit  plane  values.  Density,  p,  and  the  energy 
variable,  X,  are  available  from  integration.  Flow,  W,  is 
not  available  from  integration  and  must  be  found  using  the 
aft  boundary  condition.  At  the  exit  plane  of  the  aft 
control  volume,  static  pressure  is  specified.  Exit 
stagnation-to-static  pressure  ratio  is  tested  against 
critical  pressure  ratio  to  determine  if  the  exit  is  choked 
or  not,  as  follows: 


(89) 
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If 


^exit  Jcrit 


t 


the  exit  is  choked  and 


«in 


1.0  . 


Then, 


PS. 


in 


(Y-1)  X,„ 

TTin^  • 


(90) 


If 


^exit  J crit 


9 


the  exit  is  unchoked,  and 


Therefore, 


(91) 


Mach  number  and  static  pressure  are  now  known  for  either 
condition.  Therefore,  the  remaining  calculations  are  common 
for  either  case. 
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and 


TS 


PS. 

_ in 
in  R p.„  ' 


■\/  Y R TS.  , 
in  my'  in  ' 


(92) 

(93) 


w. 

in 


= P 


U A 

in  in  exit 


(94) 


Impulse  and  enthalpy  values  at  the  exit  and  entry 
planes  are  computed  using  Eqs.  (87)  and  (88). 

Finite  Difference  Approximation,  Central 
Difference  Form 

When  ducting  ahead  of  and  behind  the  compressor 
was  included  in  the  system  modeled,  a noncentered,  two-sided 
difference  form  of  finite  difference  approximation  was  used 
to  allow  more  accurate  representation  of  the  duct  flow 
processes  and  to  avoid  some  numerical  stability  difficulties 
experienced  in  attempts  to  apply  the  one-sided  difference 
form  to  duct-flow  cases. 

The  central  difference  algorithm  used  is  derived 
as  follows.  Let  y(z)  be  any  variable  with  continuous  third 
derivative  with  respect  to  z,  for  which  an  approximation  to 
^ IS  desired.  Referring  to  Fig.  15,  let 


and 


Az 


IP 


IP 


- z , 


(95) 


Az . 
im 


= 


im 


- *i 


(96) 
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Expanding  around  the  location, 


V = V +^l  Lz  +i 


and 


Azf  + 0 
iP 


‘^ip] 


y ■ = y • + ■5^ 
•'im  -^1  3z 


1 ^ 92^ 


Azf  + 0 
im 


[-L] 


Multiplying  Eg.  (97)  by  Eg.  (98)  by 


differencing  the  results  and  solving  for 

d Z 


gives 


where 


I2 ♦ *i  ^1  “ip  ^ip  + 0 ] - 

4^=max  . 


and  the  weighting  terms  are 


Az . 

a.  =7-^ 

im  Az 


im  [i4ip-i9i„] 


Az.  + Az. 

a = i£_ 

®i  Az. Az.‘ 

ip  im 


and 


Az . 


a . = 


im 


ip  ^^p  [az.  -Az.  1 
^ L ira  ip  J 


, (97) 


(98) 


(99) 


(100) 


(101) 
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The  sum  of  the  coefficients, 


a . + a . 

im  1 


+ a . =0 

ip 


(102) 


Multiplying  Eq.  (102)  by  y^  and  subtracting  from  Eq.  (99) 
gives 


= a . 


im 


+ a. 


iP 


(103) 


Replacing  the  spatial  derivatives  in  the  governing 
equations,  Eqs.  (2) , (4) , and  (7) , produces  the  alternate 
finite  difference  form. 


®[®]i  r„  1 f 1 

~5t l"i  - "imJ*  “ipl  "i  - “ip  ]-  ■ 110'') 

9Wi 

■5T  ' =11,  [ ] * “ip  [ 3®i  - ] 


+ F + PS 


3A 

1 5¥  . ' 
1 


(105) 


and 


[®]i 

3t 


'=iin[‘‘l-’‘i™)  •'■=ip[“i-''ta,]  flOO) 


As  was  done  with  the  one-sided  case,  the  pressure-area 
partial  derivative  in  Eq.  (105)  was  combined  with  the  axial 
force  term. 


54 


AEDC-TR-77-80 


The  solution  procedure  followed  (Fig.  16)  was  the 
same  as  used  in  the  one-sided  case  except  Eqs.  (104),  (105), 
and  (106)  were  used  instead  of  Eqs.  (63),  (64),  and  (65). 

The  bleed  flow,  forces,  shaft  wor)(  and  heat  added 
must  also  be  made  compatible  with  the  two-sided  difference 
form . Thus , 


= ^m  «®im  - ^ip  W^i  ' 

(107) 

^ ~ ^im  ^im  ^ ®ip  ^i  ' 

(108) 

"S  = - WS^  + a . p WS  . , 

(109) 

and 


Q=-a.  A.  +a.  Q. 
im  im  ip 


(110) 


Also,  on  the  forward  and  aft  boundary  planes,  an  appropriate 
one-sided  form  of  the  differential  equation  was  used.  The 
force  term  for  ducting  was  taken  as  a simple  friction  loss 
expressed  as 

F.  = - y PS^  A^  M?  . (Ill) 


The  value  of  was  adjusted  to  give  a reasonable  stag- 
nation  pressure  loss  in  the  ducting  and  to  permit 
numerically  stable  operation  of  the  model. 
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Interpretation  and  Limitations  of 
Differencing  Methods 

Differencing  algorithm  truncation  error.  The  one- 
sided algorithm  is  a first-order  approximation  (i.e.,  it 
ignores  terms  of  order  (Az)  and  above) . The  central  differ- 
ence algorithm  is  a second-order  approximation  (terms  of 
2 

oj^der  (Az)  and  above  are  ignored)  . 

Let  y again  be  any  variable  in  the  problem 
requiring  spatial  differentiation.  The  one-sided  algorithm 
is  found  from  expansion  around  a point, 

(112) 

Solving  for  the  first  partial  derivative. 


(113) 


The  truncation  error  in  the  governing  equations  is  the  error 
resulting  from  the  substitution  of  the  algorithm  for  the 
exact  form  of  the  spatial  partial  derivative.  Therefore,  an 
upper  estimate  for  the  truncation  error  for  the  one-sided 
case  is 

E = 2"  • (114) 

oZ  I. 
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The  second  partial  derivative  is  finite,  so,  as  Az  -►  0,  the 
error  vanishes  and,  referring  to  Eq.  (112) , the  approxi- 
mation is  consistent  of  order  one. 

A similar  analysis  of  truncation  error  for  the 
central  difference  algorithm  produces  an  error  estimate  of 


e 


1 

tt  — f 


(Az) 


(115) 


The  central  difference  algorithm  is  therefore  consistent  of 
order  two. 

For  the  problems  solved  using  the  nonlinear 
system  of  equations  produced  by  the  approximation,  numerical 
stability  was  investigated  by  running  various  time  step 
sizes  and  observing  the  behavior  of  the  solution  at  large 
times . 

The  system  of  equations  constituting  the  model  is 
for  t ime -dependent , one-dimensional  subsonic  flow  and  is 
hyperbolic  in  nature,  with  specified  boundary  conditions  and 
initial  conditions.  Thus,  each  integration  step  to  compute 
the  dependent  variables  at  the  new  time,  (t  + At) , must  be 
made  within  the  region  of  influence  of  the  known  values  at 
time,  t.  The  Courant-Friedrichs-Lewy  maximum  time  step 
criteria  [51] , 


Az 


At. 


min 


max 


a + U, 


max 


(116) 
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must  be  observed.  For  Az„.  = 0.25  ft,  a * 1553  ft/sec 

mrn  ' 

(value  at  TS  = 1000"R) , U = a, 

lUdX 

^^max  “ 8x10"^  sec  . 


Values  used  were  typically  one-half  of  At  . Various 

msuc 

values,  depending  on  the  nature  of  the  specific  analysis, 
were  used. 

with  consistency  and  numerical  stability,  con- 
vergence is  implied  [51] . The  best  proof  of  convergence  of 
the  numerical  solution  is  provided  by  the  comparisons  to 
experimental  results  given  in  Chapter  III. 

Disturbance  frequency  limitations.  The  differ- 
encing length  is  Az;  therefore,  no  disturbances  having 
disturbance  wave  lengths  of  order  Az  or  less  can  be 
explicitly  treated  by  the  model.  As  a working  limit, 

minimum  wave  lengths  were  restricted  to  10  Az  . There- 

zndx 

fore,  the  highest  frequency  phenomena  with  downstream 
propagation  of  pressure  disturbances  which  could  be 
explicitly  treated  by  the  model  would  be 


(U  + a) 


min 


max 


10  Az, 


max 


(117) 


where  U is  a representative  minimum  axial  velocity,  and  a is 

the  acoustic  velocity.  For  U = 300  ft/sec,  a = 1116  ft/sec, 

and  Az^ = 0.25  ft, 

max  ' 
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f = 566  Hz  . 
max 

Similar  limits  can  be  computed  for  temperature  rather  than 
pressure  oriented  disturbances. 

Steady-state  solution.  Inspection  of  the  approxi- 
mating equations,  Eqs.  (63),  (64),  and  (65),  or  Eqs.  (104), 
(105)  , and  (106) , shows  that  when  the  time  derivatives  are 
zero,  the  remaining  equations  are  the  exact  mass,  momentum 
and  energy  equations  for  steady,  one-dimensional  flow 
through  a control  volume  of  any  specified  size.  The 
governing  equations  in  the  approximate  form  could  also  have 
been  derived  from  an  integral  analysis  of  time-unsteady  flow 
through  a finite  control  volume. 
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CHAPTER  III 

EXAMPLE  PROBLEMS  WITH  COMPARISON  TO  EXPERIMENTS, 
ONE-DIMENSIONAL  MODEL 

Three  steady-state  excunples  and  three  dynamic 
examples  are  given.  Three  different  compressors,  the 
Allison  XC-1,  the  NACA-8,  and  the  General  Electric  J85-13, 
were  driven  to  their  stability  limit  for  steady  flow  to 
demonstrate  the  physical  processes  and  to  test  the  models' 
ability  to  compute  the  steady-state  stability  limit.  The 
NACA-8  compressor  was  stalled  dyneunically  by  imposed 
oscillating  inflow.  Dynamic  response  of  the  J85-13  for 
first  oscillating  inflow  and  then  for  oscillating  back 
pressure  was  computed.  Finally,  the  response  of  the  NACA-8 
compressor  to  a rapid  inlet  temperature  ramp  representing 
rocket  exhaust  gas  ingestion  is  presented. 

Three  different  compressors  were  modeled  because 
experimental  data  necessary  for  comparison  to  the  analytical 
model  did  not  exist  in  a sufficiently  complete  set  for  any 
single  compressor.  A side  benefit  of  modeling  three  com- 
pressors is  the  demonstration  of  the  generality  of  the 
analysis  method. 
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I.  STEADY-STATE  STABILITY  LIMITS 

The  Allison  XC-1  compressor  is  a four-stage 
compressor  with  design  airflow  and  stall  pressure  ratio  of 
38.7  Ibm/sec  and  4.9:1.  Details  of  the  compressor,  the 
overall  compressor  map,  and  the  stage  characteristics  used 
in  the  model  are  presented  in  Appendix  B.  The  information 
on  this  compressor  was  taken  from  [43]  . The  computational 
stations  for  the  model  are  shown  in  Fig.  17a.  No  ducting 
was  included  in  this  model,  thus,  the  one-sided  difference 
formulation  was  used. 

The  compressor  was  loaded  to  the  stability  limit 
by  increasing  the  exit  pressure  until  flow  breakdown  was 
indicated.  Inlet  stagnation  pressure  and  temperature  were 
constant.  Corrected  rotor  speed  was  90%  of  design. 

Figure  18  illustrates  the  loading  and  the  resulting  stall 
process.  Exit  static  pressure  was  ramped  at  a rapid  rate  to 
minimize  computer  time  necessary  for  the  computations.  The 
rate  was  sufficiently  slow,  however,  that  transient  effects 
were  very  small.  Airflow  at  any  stage  never  differed  from 
compressor  inlet  airflow  by  more  than  0.5%  during  loading. 

As  the  load  on  the  compressor  increased,  compressor  pressure 
ratio  increased  and  airflow  decreased  until  flow  breakdown 
was  indicated.  Figure  19  shows  some  key  compressor 
variables  on  an  expanded  time  scale  covering  the  flow  break- 
down as  computed  by  the  model.  Stage  3 flow  coefficient 
reaches  its  stall  value  at  t = 55.7  msec.  The  flow,  W^, 
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into  Stage  3 is  seen  to  begin  to  decrease  at  the  same  time. 
Shortly  afterward  the  flows  at  each  location  begin  to 
diverge  one  from  the  other  and  to  oscilla^te.  The  static 
pressure  behind  stage  3 rapidly  falls  because  the  stage's 
piimping  capacity  is  quite  limited  when  stalled.  The 
pressure  forward  of  stage  3 rapidly  increases  in  response  to 
the  sudden  blockage  to  the  flow  which  the  stalled  stage  3 
produced.  This  static  pressure  signature  is  observed 
®^P®^i***®*'tally , and  is  a common  means  employed  to  determine 
which  stage  in  a compressor  stalled  first,  causing  the 
instability  to  occur.  (See  Fig,  X-27  in  [14]  for  example.) 

Figure  20  shows  the  90%  corrected  rotor  speed 
characteristic  computed  from  the  model  compared  to  experi- 
mental measurements.  The  flow  breakdown  point  indicated  by 
the  model  agreed  well  with  the  experimental  stability  limit. 

The  model  solution  thus  indicates  the  following 
sequence  leading  to  compressor  instability.  The  compressor 
is  loaded  and  flow  reduces  until  a stage  somewhere  in  the 
compressor  reaches  its  stall  point  (its  maximum  pressure 
coefficient  point) . The  ability  of  that  stage  to  continue 
to  pump  flow  (i.e.,  to  help  support  the  pressure  gradient 
over  the  length  of  the  compressor  imposed  by  the  back 
pressure  load)  begins  to  diminish.  As  the  flow  reduction 
because  of  the  increasing  load  continues,  the  ability  of  the 
stage  to  puii^)  diminishes  even  further  because  of  the 
increased  load.  Eventually,  a point  is  reached  where  the 
flow  into  the  stage  becomes  fully  stalled,  thus,  little 
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force  on  the  fluid  is  provided  to  pump  the  flow,  and  the 
flow  path  through  the  compressor  becomes  blocked  because  of 
the  highly  separated  condition  of  the  stalled  flow  in  the 
channels.  In  terms  of  the  solution  of  the  model  equations, 
as  the  loading  is  increased,  a point  is  reached  where  no 
stable  solution  to  the  governing  equations  exists  which  will 
satisfy  the  imposed  load  (boundary  conditions);  thus,  the 
solution,  like  the  flow  through  the  compressor,  becomes 
unstable . However , the  behavior  of  the  model  quickly 
becomes  nonrepresentative  of  the  compressor  shortly  after 
flow  breakdown  because  the  stage  characteristics  are  not 
valid  for  deeply  stalled,  dynamic  situations. 

NACA-8  Compressor 

The  NACA-8  compressor  is  an  eight-stage  com- 
pressor with  design  airflow  and  stall  pressure  ratio  of 
65  Ibm/sec  and  10:1.  Compressor  details  and  stage  charac- 
teristics are  given  in  Appendix  B emd  are  taken  from  [44, 

45,  46].  The  model  included  inlet  ducting  and  combustor 
ducting,  thus,  the  central  difference  form  was  used. 
Computation  planes  are  shown  in  Fig.  17b.  Figure  21 
illustrates  the  loading  of  the  NACA-8  model  to  the  stability 
limit  by  reducing  the  turbine  nozzle  area.  The  nozzle  was 
choked.  Corrected  rotor  speed  was  90%  of  design,  inlet 
stagnation  pressure  and  temperature  were  constant.  Pressure 
ratio  increased  and  corrected  airflow  decreased  until 
instability  was  encountered.  The  process  was  repeated  at 
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70  and  807;  corrected  speed.  The  results  are  presented  in 
Fig.  22  and  are  compared  to  experimental  data.  Reasonable 
agreement  was  obtained. 

Extensive  internal  steady-state  experimental 
measurements  were  made  on  the  NACA-8  compressor  allowing 
detailed  model-to-experiment  comparisons.  Figure  23  shows 
the  variation  of  stage  flow  coefficients  near  the  stability 
limit  as  a function  of  corrected  rotor  speed.  Both  experi- 
and  model  computed  values  are  shown.  Good  agreement 
was  obtained. 

It  is  also  evident  from  Fig.  23  that  this  com- 
pressor operates  with  forward  stages  below  the  stall  value 
of  stage  flow  coefficient  at  corrected  rotor  speeds  below 
approximately  78%.  Both  experimental  and  computed  flow 
coefficient  values  for  stage  1 are  below  stall  at  70%  speed. 
This  behavior  is  in  agreement  with  previous  theories  and 
experiment  observations  (see  Chapter  I and  [21,  24] >.  An 
explanation  is  as  follows.  At  low  corrected  rotor  speeds, 
the  forward  stages  are  unable  to  deliver  the  design  pressure 
ratio.  The  flow  density  provided  to  the  aft  stages  is 
therefore  much  lower  than  that  for  which  the  compressor 
axial  area  distribution  was  designed.  Thus,  to  maintain 
steady  mass  flow  continuity,  the  flow  coefficients  (axial 
velocity)  in  the  aft  stages  are  high.  Indeed,  at  suffi- 
ciently low  corrected  rotor  speeds,  the  aft  stages  are  even 
choked.  This  limits  the  flow  through  the  compressor  to  a 
low  value.  It  also  means  that  the  aft  stage  characteristics 
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will  bo  steep  relative  to  the  forward  stages  (see  Appendix 
B) . Therefore,  at  low  corrected  rotor  speeds,  as  the  load 
on  the  compressor  (the  demanded  pressure  ratio)  is  increased, 
flow  reduces,  and  the  following  sequence  occurs.  The 
ability  of  the  aft  stages  to  increase  their  pressure  ratio 
in  response  to  the  flow  reduction  is  initially  greater  than 
the  decrease  in  pressure  ratio  occurring  in  the  forward 
stages  in  response  to  the  same  reduction  in  flow  rate.  Thus, 
the  increased  load  may  be  satisfied.  As  the  loading  con- 
tinues, however,  the  forward  stages  become  more  and  more 
stalled,  and  a point  is  eventually  reached  where  the  aft 
stages  can  no  longer  compensate.  At  that  point,  the  load 
can  no  longer  be  satisfied,  and  overall  compressor  insta- 
bility occurs. 

This  example  makes  clear  the  need  to  distinguish 
between  stage  stalling  and  overall  compressor  stalling  or 
instability.  Overall  compressor  instability  is  always 
preceded  by  stage  stall,  but  is  not  always  simultaneous  with 
stage  stall. 

General  Electric  J85-13  Compressor 

The  General  Electric  J85-13  compressor  is  an 
eight-stage  compressor  with  design  airflow  and  stall 
pressure  ratio  of  43  Ibm/sec  and  7.7:1.  The  compressor 
employs  a set  of  variable  inlet  guide  vanes  ahead  of  the 
first-stage  and  compressor  bleed  from  the  second-,  third-, 
and  fourth-stage  stator  cases.  The  guide  vanes  and  bleed 
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valve  are  scheduled  by  the  engine  control  as  functions  of 
corrected  rotor  speed.  (Additional  bias  for  inlet  tempera- 
ture is  also  sometimes  provided.)  The  schedules  used  in  the 
model  of  this  study  are  given  in  Appendix  B.  Stage  charac- 
teristics and  additional  compressor  details  are  also  given 
in  Appendix  B and  are  taken  from  [47] . 

The  steady-state  stability  limit  and  corrected 
speed  characteristics  were  computed  for  80  and  87%  corrected 
rotor  speed  and  are  compared  to  experimental  results 
(Fig.  24).  Agreement  was  good.  The  compressor  model  was 
loaded  by  reducing  turbine  nozzle  area.  The  nozzle  was 
choked . The  model  included  ducting , thus  the  central 
difference  formulation  was  used. 

II.  DYNAMIC  STABILITY  LIMITS 

Oscillating  Inflow  Loading,  NACA-8  Compressor 

Figure  25a  shows  the  oscillation  of  the  inlet 
total  pressure  imposed  on  the  NACA-8  compressor  model. 

Inlet  total  temperature  and  turbine  nozzle  area  were 
constant.  Corrected  rotor  speed  was  70%.  The  pressure  was 
oscillated  at  20  Hz  and  the  amplitude  was  slowly  increased 
until  compressor  instability  occurred.  Overall  compressor 
pressure  ratio  responded  as  shown  in  Fig.  25b.  For  con- 
venience, the  model  was  started  at  90%  speed  (initial 
conditions  were  previously  calculated  for  90%  speed) , 
rapidly  decelerated  to  70%  speed,  stabilized,  and  then  the 
oscillation  was  imposed  at  t = 0.15  second. 
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The  pressure  ratio  oscillations  grew  with  the 
imposed  pressure  oscillations  until  instability  was 
encountered.  The  response  of  the  airflow  into  the  com- 
pressor is  shown  in  Fig,  25c.  Initially,  the  flow  variation 
is  near  sinusoidal,  but,  as  the  stages  are  driven  into  and 
out  of  their  stalled  regions,  indications  of  flow  breakdown 
begin  to  appear.  Finally,  the  flow  completely  diverges  to 
zero.  The  phase  shifts  and  amplitude  variations  of  the  flow 
properties  through  the  compressor  due  to  the  forced  oscil- 
lations cause  the  stage  flow  coefficients  to  oscillate. 
Figures  25d  through  25g  show  the  flow  coefficient  dynamic 
variations.  Also  indicated  are  the  individual  stage  stall 
values  and  the  value  of  each  flow  coefficient  at  which  the 
compressor  steady-state  stability  limit  was  reached  at  70% 
rotor  speed.  Excursions  well  below  both  values  occur 
dynamically  before  total  flow  breakdown  is  experienced. 

This  behavior  is  in  qualitative  agreement  with  experimental 
observations  reported  in  [5,  37] . Examinations  of  the  flow 
coefficient  variations  indicate  that  stage  7 was  the  stage 
which  eventually  became  totally  unstable  and  caused  com- 
pressor instability  (Fig.  25g) . 

Interpretation  of  the  model  results  suggests  the 
following  physical  process  for  instability  with  oscillating 
inflow.  The  stage  flow  coefficients  are  alternately  driven 
into  and  out  of  their  stalled  regions.  The  stages  must  be 
held  in  the  adversely  loaded  condition  for  a sufficiently 
long  period  of  time  for  the  flow  breakdown  to  occur.  Thus, 
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at  a given  frequency,  the  amount  of  time  a stage  (or  group 
of  stages)  will  be  in  the  highly  loaded  condition  depends  on 
the  amplitude  of  the  imposed  oscillation.  Once  a stage  (or 
group  of  stages)  is  transiently  held  in  the  adversely  loaded 
condition  for  a sufficient  time  period,  the  breakdown 
process  is  locally  like  that  which  occurs  in  steady-state 
loading . 

The  extent  to  which  the  flow  coefficients  oscil- 
late depends  on  the  dynamic  response  (phase  angle  and 
amplitude  ratio) , which  depends  on  oscillation  frequency. 

The  loading  process  of  Fig.  25  was  repeated  for  frequencies 
up  to  320  Hz.  The  amplitude  of  imposed  oscillation  required 
to  drive  the  compressor  to  the  stability  limit  is  shown  in 
Fig.  26.  As  frequency  increases,  the  amplitude  required  for 
instability  decreases.  This  occurs  because,  as  frequency 
increases,  the  phase  lags  and  amplitude  attenuations  across 
the  stages  increase;  thus,  the  flow  coefficient  oscillations 
are  larger  for  a given  imposed  oscillation  amplitude  at 
higher  frequencies.  (A  fundamental  discussion  of  the  phase 
shift  and  amplitude  ratio  effect  is  given  in  [24].) 

A local  minimum  occurs  in  Fig.  26  near  80  Hz. 

This  occurs  because,  below  80  Hz,  aft  stages  in  the  com- 
pressor are  critical,  but,  at  frequencies  higher  than  80  Hz, 
the  amplitude  attenuation  across  the  forward  stages 
increases  to  the  extent  that  it  protects  the  aft  stages,  and 
a new  stage,  farther  forward  in  the  compressor,  becomes 
critical . 
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The  computations  of  Fig.  26  were  made  both  with 
and  without  the  unsteady  cascade  correction  to  the  stage 
characteristics.  First-stage  reduced  frequency  is  indicated 
at  each  value  of  oscillation  frequency.  Only  small 
influences  were  experienced  for  this  compressor  in  oscil- 
latory inflow  because  of  the  relatively  low  reduced 
frequencies.  Compressors  and  fans  with  larger  blades  would 
experience  the  influence  of  the  correction  at  somewhat  lower 
frequencies. 

No  experimental  data  were  available  for  this  com- 
pressor for  direct  comparison  of  stability  limits.  The 
general  trend  with  frequency  and  the  amplitude  magnitudes 
are  in  agreement  with  the  experiment  of  [37],  however. 

Oscillating  Inflow  Dynamic  Response, 

J85-13  Compressor 

As  indicated  in  the  previous  section,  the  ability 
of  the  model  to  properly  compute  stability  limits  in  oscil- 
latory flow  depends  on  its  ability  to  compute  the  dynamic 
response  of  the  various  stages  in  the  compressor  to  the 
imposed  oscillations.  Two  experiments  were  carried  out  by 
NASA  using  the  J85-13  turbojet  engine  to  measure  the 
un stalled  dynamic  response  of  the  compressor  to  imposed 
inflow  oscillations  [52,  53].  Figure  27  is  a sketch  of  the 
test  arrangements . One  test  was  conducted  in  a supersonic 
wind  tunnel  with  an  axisymmetric  mixed  compression  inlet 
installed  (Fig.  27a) . Engine  entry  total  pressure 
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oscillations  were  produced  by  oscillating  the  inlet  bypass 
doors,  causing  the  terminal  normal  shock  in  the  inlet  to 
oscillate,  thus  producing  the  downstream  total  pressure 
oscillations.  The  other  test  arrangement  (Fig.  27b)  used  an 
system  in  a section  of  ducting  ahead  of  the  engine. 
The  air jet  system  operates  by  blowing  high  velocity  air 
counter  to  the  primary  airstream.  Planar  oscillations  in 
engine  entry  pressure  were  produced  by  oscillating  the  air- 
jet  counterflow.  In  both  cases,  the  oscillation  amplitude 
was  approximately  10%,  peak-to-peak. 

The  J85-13  compressor  model  was  subjected  to 
imposed  total  pressure  oscillations  at  the  forward  control 
volume  plane,  thus  simulating  the  inlet  experiment  of 
Fig.  27a.  The  test  setup  of  Fig.  27b  differs  from  the  model 
and  the  inlet  experiment,  in  that  the  forward  isolation 
plane  provided  by  the  supersonic  flow  in  the  throat  of  the 
inlet  is  not  in  the  air jet  test  setup.  The  inlet  and  model 
dynamic  systems  consisted  of  the  ducting  downstream  of  the 
normal  shock  wave  to  the  choked  turbine  nozzle  in  the  engine. 
The  airjet  dynamic  system  includes  all  of  the  test  cell 
ducting  forward  of  the  engine  (until  a choked  valve, 
venturi,  or  other  forward  isolation  point  is  reached)  down 
to  the  choked  turbine  nozzle  in  the  engine. 

Figure  28  shows  the  variation  of  pressure  amplitude 
ratio  and  phase  angle  for  a 40-Hz  oscillation  through  the 
compressor.  Model  computations  are  compared  to  the  NASA 
data.  The  calculated  trends  are  in  good  agreement  with 
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experiment.  Normalized  magnitude  differs  by  approximately 
10%  and  phase  angle  by  five  degrees.  The  agreement, 
although  reasonable  for  dynamic  response  data,  could  be 
improved  by  closer  matching  of  engine  operating  condition 
between  model  and  experiment.  Also,  the  J85-13  engine  run 
in  the  wind  tunnel  was  likely  not  the  same  engine  from  which 
the  stage  characteristics  were  taken;  thus,  engine-to-engine 
variations  are  present.  Further,  the  measurement  uncertain- 
ties present  in  the  experimental  data  (quoted  as  ±3%  on 
pressure  measurement  and  ±8  degrees  on  phase  angle  in  [52] ) 
may  affect  the  comparison. 

More  detailed  dynamic  response  results  are  shown 
in  Fig.  29.  Figures  29a  and  29b  represent  the  frequency 
response  across  one  stage;  Figs.  29c  and  39d  across  approxi- 
mately half  of  the  compressor;  and  Figs.  29e  and  29f  across 
the  entire  compressor.  Trend  and  level  agreement  between 
the  model  computations  and  the  inlet  results  were  reasonable, 
but,  once  again,  could  be  improved  by  closer  matching  of  the 
actual  experimental  hardware.  In  particular,  the  compressor 
discharge  agreement  suffered  because  the  combustor  volume  of 
the  actual  engine  was  larger  than  that  modeled.  The  com- 
bustor in  the  model  was  shortened  during  steady-state 
stability  studies  to  keep  the  number  of  computation  planes 
to  a minimum,  and  was  not  increased  for  the  dynamic  response 
study. 

Comparisons  of  the  model  and  engine-with-inlet 
results  with  the  airjet  results  show  considerable 
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differences,  especially  near  40  Hz.  The  most  likely  cause 
for  the  difference  is  the  fact  that  the  two  dynamic  systems 
are  not  equivalent. 

The  model  configuration  used  in  the  computations 
was  not  tailored  in  detail  to  the  experimental  configuration 
tested.  The  comparison  shows  that  to  obtain  close  agreement 
dynamically,  close  attention  to  matching  configuration  must 
be  observed.  The  comparisons  also  show  that  considerable 
care  must  be  exercised  in  dynamic  experimental  work  to 
assure  dynamic  equivalency  of  systems  (e.g.,  same  isolation 
points,  volumes  and  geometry) . 

The  results  also  suggest  a potentially  powerful 
coupling  of  experimental  work  and  analytical  modeling.  For 
example,  if  it  were  not  possible  (or  desired)  to  physically 
dynamically  isolate  the  airjet  setup  in  the  test  cell  to 
simulate  a given  inlet-engine  combination,  the  test  setup, 
including  the  airjet  system  and  the  compressor,  could  be 
modeled  and  the  model  refined  using  experimental  results  to 
closely  match  the  test  setup  measurements.  Then,  the 
analytical  model  inlet  configuration  could  be  modified  to 
match  any  desired  configuration,  and  the  desired  information 
derived  from  the  model. 

Oscillating  Compressor  Discharge  Pressure, 

J85-13  Compressor 

Fuel  flow  to  the  J85-13  engine  run  with  the  super- 
sonic inlet  (Fig.  27a)  was  oscillated  to  determine  the 


72 


AEDC-TR-77-80 


propulsion  system's  dynamic  response  to  oscillations  induced 
within  the  engine.  The  oscillations  produced  compressor 
discharge  oscillations  of  approximately  10%.  The  inlet  was 
assumed  to  remain  started,  thus,  the  model  dynamic  system 
boundaries  were  the  same  as  for  the  oscillating  inlet  bypass 
door  case.  Computations  were  made  using  the  J85-13  model  by 
adding  heat  to  the  combustor  control  volumes  in  an  oscil- 
latory manner  to  produce  peak-to-peak  oscillations  in 
compressor  back  pressure  of  approximately  10%. 

Figure  30  provides  a comparison  between  computed 
and  experimental  amplitude  ratio  and  phase  angle  variation 
through  the  compressor  at  a frequency  of  40  Hz.  Trend  and 
level  agreement  were  reasonable.  Because  of  the  high  com- 
pressor rotor  inertia,  at  40  Hz,  rotor  speed  cannot  respond 
to  the  imposed  fuel  flow  oscillations;  thus,  the  effects  of 
the  oscillation  are  from  aerodynamic  upstream  propagation 
only.  The  aft-stage  characteristics  (stages  4 through  8) 
(Appendix  B)  are  quite  steep,  and  at  this  condition,  the 
stages  are  operating  on  the  steep,  negative  slope  part  of 
the  curves.  Therefore,  perturbations  in  the  aft-stage 
airflow  (or  flow  coefficient)  will  cause  large  oscillations 
in  stage  pressure  ratio.  Figure  30  indicates  the  magnifi- 
cation of  the  induced  oscillation  at  the  stage  4 entry. 
However,  stages  1,  2,  and  3 are  operating  in  the  flow  coef- 
ficient range  of  0.5  to  0.55  (i.e. , on  the  shallow,  positive 
slope  part  of  the  curves  for  stages  2 and  3,  and  the 
shallow,  negative  slope  range  for  stage  1) . As  a result. 
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the  forward  stages  oscillate  pressure  ratio  very  little  in 
response  to  a given  airflow  oscillation.  Therefore,  the 
forward  stages  attenuate  the  upstream  propagating  dis- 
turbance to  the  point  that  it  is  barely  measurable  at  the 
compressor  inlet  (Fig.  30a) . Thus,  it  becomes  clear  that 
the  stage  characteristic  shapes  and  match  points  play  a very 
important  role,  not  only  in  the  steady-state  stability  of 
the  compressor,  but  also  in  the  way  disturbances  are  propa- 
gated, thus,  in  the  transient  and  dynamic  stability 
characteristics  of  the  compressor. 

Figures  31a  and  31b  present  the  dyneunic  response 
to  oscillating  fuel  flow  across  the  compressor.  Experi- 
mental results  are  presented  along  with  three  different  sets 
of  model  computations.  The  first  set  is  for  constant  rotor 
speed  and  constant  and  T^^.  The  second  set  is  for 
variable  rotor  speed  and  constant  and  T^^.  The  third  set 
is  for  variable  rotor  speed,  variable  Pj^  and  constant  T^. 

The  three  sets  were  computed  to  account  for  the  three 
dynamic  effects  caused  by  oscillating  fuel  flow.  At  fre- 
quencies above  approximately  20  Hz,  the  rotor  speed  cannot 
respond  to  the  oscillating  fuel  flow;  thus,  as  already 
discussed,  upstream  aerodynamic  propagation  is  the  chief 
mechanism  present.  At  frequencies  below  20  Hz,  rotor  speed 
can  respond  to  the  oscillations  in  fuel  flow,  thus,  the 
stage  flow  coefficients  must  now  vary  mechanically  also 
(i.e.,  wheel  speed  oscillates).  Further,  because  wheel 
speed  varies,  compressor  airflow  demand  also  must  vary. 
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Inlet  airflow  varies  as  a result  and  causes  the  normal  shock 
wave  to  oscillate.  This  causes  the  inlet  total  pressure, 
to  oscillate. 

Engine  rotor  speed  response  was  represented  in  the 
model  by  a simple  first-order  differential  equation  simu- 
lating the  engine  control  and  rotor  dynamic  response.  An 
oscillation  in  at  the  control  volume  forward  boundary, 
was  imposed  to  simulate  the  shock  wave  motion.  As  seen  in 
Fig.  31,  reasonable  agreement  with  experiments  was  obtained 
only  when  all  three  of  these  effects  were  included.  Because 
of  the  complicated  nature  of  the  three  disturbance  paths, 
agreement  between  model  computations  and  experimental 
results  at  internal  compressor  locations  was  poor.  Agree- 
ment could  be  improved  by  more  exact  modeling  of  the 
specific  engine  under  consideration. 

This  example  makes  clear  the  importance  of 
auxiliary  perturbation  propagation  paths  which  must  be  con- 
sidered, and  which  can  have  a strong  influence  on  dynamic 
response,  thus  on  stability. 

Rapid  Inlet  Temperature  Rcunp  to  Stall, 

NACA-8  Compressor 

The  NACA-8  compressor  model  was  subjected  to  a 
rapid  rise  in  inlet  temperature  simulating  ingestion  of  hot 
gases  from  gun  or  rocket  fire.  Mechanical  rotor  speed,  N, 
and  were  held  constant.  The  initial  corrected  rotor 
speed  was  70%.  Figure  32a  shows  the  imposed  temperature 
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ramp.  The  engine  entry  temperature  responded  as  shown  in 
the  figure.  Approximately  0.01  second  was  required  for  the 
temperature  pulse  to  be  convected  to  the  engine  face. 

Figure  32b  shows  the  rapid  reduction  in  corrected  rotor 
speed  caused  by  the  temperature  ramp.  Figure  32c  shows  the 
stage  flow  coefficient  variation  caused  by  the  temperature 
ramp.  The  forward  stages  are  rapidly  driven  downward  by  the 
reduction  in  airflow  accompanying  the  rapid  temperature  rise. 
The  resulting  increase  in  forward  stage  loading  immediately 
prior  to  stall  initially  holds  the  mid-stage  reduction  to  a 
minimum,  and  actually  unloads  the  aft  stages.  A point  is 
reached,  however,  where  the  forward  stages  stall  and  can  no 
longer  support  the  load  imposed  on  them.  Flow  breakdown 
then  occurs  and  is  propagated  back  through  the  compressor. 

In  the  example  shown,  the  compressor  may  have 
recovered  from  the  instability  with  the  removal  of  the 
temperature  pulse.  On  an  actual  engine,  combustor  flameout 
would  probably  have  occurred  as  a result  of  the  compressor 
instability  which  was  initiated. 

A series  of  ramps  of  differing  ramp  height  was 
computed  at  a ramp  rate  of  3600“R/second.  The  results  are 
compared  in  Fig.  33  to  a scatter  band  for  three  different 
engines  from  [16]  . Although  Fig.  33  is  far  from  an  exact 
comparison,  the  results  do  indicate  proper  trend  and 
magnitude . 
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Gas  composition  was  considered  constant  for  these 
calculations.  In  the  real  experimental  case,  additional 
effects  may  be  present  because  of  the  gas  composition 
changes  associated  with  ingestion  of  combustion  products. 
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CHAPTER  IV 

DEVELOPMENT  OF  THREE-DIMENSIONAL  MODEL 
FOR  DISTORTED  FLOWS 

The  concepts  used  to  develop  the  one-dimensional 
model  can  be  extended  to  allow  computation  of  spatial  dis- 
tortion effects  on  stability.  Figure  34  depicts  the  control 
volume  arrangement  used  for  the  distortion  model.  As  was 
noted  in  Chapter  II,  the  one-sided  finite  difference  form  of 
the  governing  equations  could  have  been  derived  from  a 
finite  control  volume  integral  analysis  rather  than  a 
differential  analysis.  For  greater  simplicity  in  developing 
the  distortion  model,  the  finite  control  volume,  integral 
approach  is  used. 

The  development  of  the  distortion  model  is 
described  as  follows.  The  governing  equations  are  written. 
Simplifications  to  allow  distortion  effects  to  be  computed 
with  a minimum  amount  of  internal  radial  and  circumferential 
flow  detail  are  described.  Special  physical  considerations 
for  radial  distortion  effects  and  for  circumferential  dis- 
tortion effects  are  presented.  The  method  of  solution, 
closely  paralleling  that  used  for  the  one-dimensional  model 
is  then  described. 
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I.  GOVERNING  EQUATIONS 

Figure  34  shows  a finite  control  volume  covering 
an  arc  sector  of  a concentric  ring  of  the  compressor  annulus 
extending  axially  into  the  compressor  one  stage  length.  The 
compressor  is  divided  into  control  volumes  as  indicated  in 
Fig.  34.  The  subscript  notation  used  to  identify  the 
location  of  a control  volume  is  also  shown.  Subscript  "i" 
indicates  an  axial  position,  "j"  a radial  position,  and  "k" 
a circumferential  position.  Control  volume  i,j,k  is  there- 
fore the  control  volume  whose  axial-facing  entry  plane  is 
on  plane  i,  and  whose  axial-facing  exit  plane  is  on  plane 
ip  = i + 1.  Its  lower  boundary  surface  approximately 
facing  the  radial  direction  is  on  surface  j,  its  upper  face 
is  on  surface  jp  = j + 1.  Similarly,  the  control  volume 
surface  normal  to  the  circumferential  direction  is  on 
surface  k and,  advancing  counterclockwise,  its  other  surface 
normal  to  the  circumferential  direction  is  on  plane 
kp  = k + 1.  The  compressor  inlet  is  at  i = 1 , the  exit  at 
i = in.  The  annulus  inner  surface  is  at  j = 1,  the  outer 
surface  is  at  j = jn.  Circumferentially,  k = 1 is  the  same 
plane  as  k = kn  + 1.  The  compressor  is  therefore  divided 
into  inm  • jnra  • kn  control  volumes  where  inm  = (in-1)  and 
jnm  = (jn-1) . For  the  example  of  Fig.  34,  the  compressor  is 
divided  into  72  control  volumes. 

The  compressor  is  divided  radially  on  an  equal 
axially-facing  area  basis.  That  is, 
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total . 

AZ  . ..  = , 

ijk  3nm  • kn  ' 


(118) 


where  is  the  total  axial  facing  area  at  any  plane, 

■‘‘i 

i.  Thus,  all  control  volumes  having  the  scune  index,  i,  are 
of  equal  volume  because  their  lengths  are  equal. 

Figure  35  illustrates  a single  control  volume. 

For  the  purpose  of  this  analysis,  the  control  volume  is 
assumed  approximately  rectangular  on  each  face  and  therefore 
is  a parallelepiped  with  the  coordinate  system;  z,  the  axial 
direction,  r,  the  radial  direction,  and  s,  a coordinate 
along  the  arc  length  of  the  control  volume.  The  governing 
equations  for  each  control  volume  may  be  written  as  follows. 


Mass 


WZ . . , + WR . . , + WC . . , = 
13k  13k  13k 


mass  flow  entering  the  control 
volume  per  unit  time 


WZ . . . + WR . . , + WC . . , 

ip3k  i3pk  i]kp 


mass  flow  leaving  the  control 
volume  per  unit  time 


cd (vol) 


(119) 


time  rate  of  increase  of  mass 
"stored"  in  the  volume 
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where , 


WZ.  ..  , WR.  ..  , WC.  .. 

13k  ' 13k  ' 13k 


are  the  mass  flows  into  the  axially-facing,  radially-facing, 
and  circumferentially— facing  surfaces,  respectively. 

Axial  Momentum 


* “1  jk  ■ ”ipjk  “ipjk  1 


axial  forces  acting  on  control  volume 


”^ipjk  ”ipjk  *^^i  jpk  ^i jpk  ^^ijkp  ^ijkp 


axial  momentum  leaving  control  volume 
per  unit  time 


- W='i3kUi3k-«®i3k“i3k-«^ijk^i3k 


axial  momentum  entering  control  volume 
per  unit  time 


i 


^ I (pU)d(vol) 

fvol.  .. 

13k 


y (120) 


time  rate  of  increase  of  axial  momentum 
"stored"  in  the  control  volume 
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The  term,  FZ,  is  the  force  of  the  blading  and 
cases  acting  on  the  fluid  in  the  axial  direction.  In 
Eq.  (120) , a projected  area  convention  is  used.  The  term, 
ARZ,  means  the  projected  area  in  the  axial  (Z)  direction  of 
the  radially— facing  (R)  surface.  Similarly,  the  projected 
area  in  the  axial  direction  of  the  circumferentially-facing 
area  is  ACZ,  and  is  zero  because  of  the  way  the  control 
volume  was  defined. 

The  bar  above  certain  quantities  of  Eq.  (120) 
indicates  that  these  terms  actually  represent  surface 
integrals.  For  example. 


WS 


ARZ. 


ijk  ""“ijk 


PS  d(ARZ)  . 
ijk 


(121) ' 


That  is,  the  barred  quantity  indicates  an  integration  of 
pressure  over  the  projected  area  of  interest.  Similarly, 


WR 


ijpk 


pvU  d(ARR)  . 
ijpk 


(122) 


Here,  the  barred  quantity  indicates  integration  over  the 
radial-facing  surface  area  of  the  axial  momentum  being 
carried  across  the  surface  by  the  radial  flow.  Implicit  in 
Eq.  (120)  is  the  assumption  that  the  flow  properties  are  all 
uniform  over  the  axially-facing  surfaces  of  any  given 
control  volume  because  no  surface  integration  is  indicated 
over  those  surfaces. 
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Radial  Momentum 


FR.  ..  + PS  . ..  ACR.  + PS.  ACR. 
i]k  ijk  ijk  ijkp  i3kp 


+ PS  . . , ARR ...  - PS  . . , ARR . . , 
i]k  i]k  ijpk  iDPk 


radial  forces  acting  on  control  volume 


WZ.  .^v.  ..+WR..,v..,+WC...  w... 

ipjk  ipjk  i]pk  i3pk  ijkp  i^kp 


radial  momentum  leaving  control  volume 


WZ . . . V . . , + WR . . . V . . , + WC . . . w . . , > 

13k  ijk  13k  ijk  13k  13k 


radial  momentum  entering  control  volume 


( pv) d (vol) 


time  rate  of  increase  of  radial  momentum 
"stored"  in  the  control  volume 


(123) 


The  force,  FR,  includes  centrifugal  forces  and  viscous 
forces  acting  radially  on  the  fluid  in  the  control  volume. 
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Circumferential  Momentum 


FC  . . , + . . . ACC  . . . - PS  . ..  ACC. 

i]k  ijk  ijk  ijkp  ijkp 


circvunferential  forces  acting  on 
control  volume 


WZ  . . . w . . , + WR . . , w . . , + WC ...  w . . , 
ip]k  xpjk  iDpk  iDpk  i]kp  i]kp 


circumferential  momentum  leaving  control 
volume  per  unit  time 


WZ.  W...-WR...  w..,  -WC...  w. 
i]k  i]k  i]k  i]k  idJ«  ijk 


^ . (124) 


circumferential  momentum  entering  control 
volume  per  unit  time 


(ow)d(vol) 


vol . 
i:k 


time  rate  of  increase  of  circumferential 
mcmentum  "stored"  in  the  control  volume 


The  force,  FC , includes  Coriolis  and  viscous  forces  acting 
on  the  fluid  in  the  control  volume. 
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Energy 


HZ.  ..  + HR.  .,  + HC.  .. 
13k  13k  13k 


enthalpy  entering  the  control 
volume  per  unit  time 


+ WS  + Q 


shaft  work  done  and  heat  added 
to  fluid  in  control  volume  per 
unit  time 


(125) 


enthalpy  leaving  the  control 
volume  per  unit  time 


e + 


2 2 2 
U + V + w 


vol 

13k 


d (vol) 


time  rate  of  increase  of  energy 
"stored"  in  the  control  volume 


In  Eq.  (125) , the  barred  quantities  represent  integrations 
over  surfaces  similar  to  those  of  the  momentum  equations. 
For  example, 


pv  d (AR) 

. 


(126) 
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II.  SIMPLIFICATIONS 

In  principle,  the  governing  equations,  Eqs.  (119) , 
(120)  , (123)  , (124)  , and  (125)  , could  be  solved  directly  if 
all  terms  entering  the  equations  could  be  evaluated.  There 
are  five  differential  equations  for  each  control  volume; 
therefore,  for  the  excunple  of  Fig.  34,  simultaneous  solution 
of  360  differential  equations  would  be  required.  This,  in 
itself  is  not  prohibitive,  but,  as  finer  divisions  of  the 
compressor,  or,  as  larger  compressors  are  considered,  the 
computer  storage  and  time  requirements  make  the  analysis 
inconvenient  to  use.  Also,  direct  solution  of  the  radial 
and  circumferential  momentvun  equations  requires  knowledge  of 
equivalent  radial  and  circumferential  body  forces  not 
readily  determinable.  Further,  numerical  stability  becomes 
more  and  more  of  a problem  as  the  level  of  complexity  and 
number  of  equations  are  increased. 

The  aim  of  this  analysis  is  to  compute  the 
influence  of  distorted  flows  on  compressor  stability.  This 
goal  can  be  achieved  with  simplifications  to  the  governing 
equations  which  circumvent,  at  least  to  an  extent,  the 
problems  listed  above  and  produce  a practical  model  for  dis- 
tortion analysis.  The  simplifications  are  as  follows. 

Radial  and  Circumferential  Flows 

Each  control  volume  of  the  compressor  is  assumed 
to  be  representable  as  a compressor  element  which  is  in 
radial  and  circumferential  equilibrium  with  adjacent  control 
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volumes  when  the  compressor  is  running  with  uniform  flow  at 
its  entry  plane  and  uniform  throttling  at  the  exit.  This 
means  radial  and  circumferential  crossflows  are  zero  with  no 
distortion  present.  Further,  it  is  assumed  that  with  dis- 
tortion present,  the  radial  and  circumferential  flow 
velocities  are  no  greater  than  an  order  of  magnitude  less 
than  the  axial  velocity. 

The  circumferential  and  radial  crossflows  into  the 
control  volume  were  computed  using  empirical  relationships 
in  place  of  the  circumferential  and  radial  momentum 
equations.  Algebraic  expressions,  relating  the  crossflows 
to  the  static  pressure  differences  between  volumes,  were 
used.  This  substitution  reduced  the  number  of  differential 
equations  requiring  solution  from  five  times  to  three  times 
the  number  of  control  volumes. 

Circumferential  flows.  Figure  36  illustrates  a 
pair  of  control  volumes  circumferentially  adjacent  but  with 
different  static  pressures.  Looking  down  the  blades  in  a 
spanwise  direction,  the  flow  path  for  the  circumferential 
flow  IS  through  the  blade  row  gaps.  The  blade  row  gaps,  in 
effect,  form  a series  of  orifice-like  restrictions  to  the 
circumferential  flow.  The  circumferential  flow  is  there- 
fore computed  by  assuming  that  the  static  pressure  in  the 
high  pressure  region,  is  a "reservoir"  pressure 

driving  flow  to  a lower  pressure,  across  a series  of 
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orifices.  The  sequence  of  equations  comprising  the  approxi- 
mating algorithm  is  as  follows . 


(127) 


(128) 


WC 


gap 


D,  V A 

^low  pressure  gap  gap 

region  effective 


(129) 


The  approximations  of  using  low  pressure  region  acoustic 
velocity  and  density  are  justified  because  the  actual 
differences  between  adjacent  control  volumes  are  not  large. 
The  effective  gap  area  is  given  by 


*gap 

effective 


CXFC 


A 

gap 

physical 


(130) 


where  CXFC  is  the  circumferential  crossflow  coefficient. 

For  most  cases,  a nominal  value  of  CXFC  = 0.6  was  used.  The 
best  way  to  determine  the  coefficient  value  would  be  to 
empirically  calibrate  it  using  compressor  test  data  to  pro- 
vide observed  crossflow  values  with  circumferential 
distortion  applied  to  the  compressor.  The  physical  gap  area 
used  included  the  space  between  the  rotor  and  the  stator 
internal  to  each  control  volume  plus  the  gap  between  the 
trailing  edge  of  the  stator  upstream  of  the  control  volume 
and  the  rotor  of  the  control  volume.  For  the  first  stage. 
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the  gap  area  was  doubled  to  provide  an  accounting  for  the 
additional  circumferential  flow  adjustment  which  occurs 
upstrecun  of  the  compressor  first  rotor.  Experimental  data 
is  given  in  [54]  which  indicates  the  extent  of  the  cir- 
cumferential flow  redistribution  occurring  at  the  compressor 
entry.  The  process  is  described  in  detail  in  [38] . The 
impact  of  crossflow  magnitude  on  the  solution  is  indicated 
in  Chapter  V. 

The  circumferential  velocity  component  across  the 
face  of  a control  volume  was  computed  as 


13k  2 


(131) 


The  averaging  is  performed  because  the  "gap"  values  apply 
at  the  circumferential  boundaries  of  the  control  volume. 

Radial  flows.  Figure  37  shows  a pair  of  radially 
adjacent  control  volumes  with  different  static  pressure. 
Viewed  looking  axially  into  the  compressor,  the  blades  form 
channels  in  which  radial  flow  may  occur.  In  the  actual 
compressor,  strong  secondary  flow  fields  caused  by  rotation 
and  viscous  effects  exist  [55]  . A considerable  amount  of 
mass  is  carried  back  and  forth  across  the  control  volume 
radial  boundary  on  a scale  usually  smaller  than,  or  on  the 
order  of,  the  control  volume  size.  For  the  purposes  of  this 
analysis,  it  is  the  net  flow  exchange,  caused  by  distortion. 
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which  is  desired.  An  approximate,  empirical  form  of  the 
radial  momentum  equation  is  used. 


WR, 


net 


= exFR 


PS 


high 


- PS 


low 


+ FR 


(132) 


The  form  of  the  equation  is  a simplification  of  the  full 
radial  momentum  equation  produced  by  neglecting  the  time- 
dependent  term  and  flux  terms  across  all  but  the  radial- 
facing  control  volume  boundaries.  The  force  term,  FR,  is 
representative  of  the  sum  of  viscous  and  centrifugal  forces. 


FR  = F . + p vol  u r 

viscous 


(133) 


Because  WR^^^  is  the  radial  crossflow  caused  by  distortion, 

the  force,  FR,  is  adjusted  so  as  to  cause  WR  . to  be  zero 

net: 

with  the  compressor  running  with  no  distortion.  This  is,  in 
a sense,  equivalent  to  imposing  the  radial  equilibrium 
condition. 


= 

3r 


2 

p w r , 


(134) 


often  used  in  three-dimensional  analysis  of  compressor 
flows  [20]  . (In  this  application,  however,  FR  conceptually 
includes  viscous  forces  as  well  as  centrifugal  forces.)  If 
the  stage  characteristics  of  each  radially  adjacent  control 
volume  were  the  same,  with  no  distortion  present,  FR  would 
be  zero.  If  radial  variation  of  stage  characteristics  is 
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considered  (treated  later  in  this  chapter) , FR  is  set  to 
cause  to  be  zero  with  no  distortion  present. 

The  term,  CXFR,  is  the  radial  crossflow  coef- 
ficient, and  was  determined  by  subjecting  the  model  to 
radial  inlet  distortion  and  adjusting  the  coefficient  to 
permit  the  stage  exit  radial  static  pressure  gradient  to 
vary  as  much  as  2%  from  the  undistorted  value.  This  pro- 
cedure was  adopted  because  experimental  data  on  radial 
distortion  indicate  that  although  the  induced  static 
pressure  gradients  relax  as  the  flow  progresses  through  the 
compressor,  a slight  gradient  does  remain  [561.  The  upper 
limit  of  CXFR  was  set  by  model  numerical  stability  con- 
siderations. A value  of  CXFR  = 0.01  was  used  for  most  cases 
computed  in  this  study.  For  modeling  a specific  compressor, 
the  value  of  CXFR  would  best  be  determined  through  tests  of 
a compressor  with  radial  distortion  imposed  and,  adjustment 
of  the  coefficient  made  to  match  the  observed  change  in 
stage  exit  flow  distributions  caused  by  the  radial 
distortion . 

The  radial  velocity  across  the  boundary  is 
computed  from. 


VR 


WR 


net 


o A " 

^low  channel 


(135) 
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The  radial  velocity  component  across  the  face  of  the  control 
volume  was  computed  as 


^ijk 


VR.  ..  + VR.  . , 
i]k  ' i]pk 


(136) 


III.  SPECIAL  CONSIDERATIONS  FOR 
DISTORTION  EFFECTS 

Radial  Distortion 

Figure  38  illustrates  the  key  physical  processes 
occurring  when  a radially  distorted  flow  approaches  a rotor 
blade.  On  the  left  of  Fig.  38,  with  uniform  flow,  con- 
sidering the  "design  A"  case,  the  tip  and  hub  angles  of 
attack  and  stage  characteristic  operating  points  are  stable, 
near  design  conditions.  On  the  right  side  of  Fig.  38,  a hub 
radial  distortion  (low  velocity  or  total  pressure  in  the  hub 
region)  case  is  shown.  For  case  "A,"  the  tip  section 
operates  at  a lower  angle  of  attack  and  higher  flow  coef- 
ficient while  the  hub  becomes  highly  loaded,  and,  in  case 
"A,"  is  even  stalled. 

In  the  "design  "B"  case,  the  tip  section  in  uni- 
form flow  is  highly  loaded  while  the  hub  section  is  only 
lightly  loaded.  Introduction  of  the  radial  distortion 
pattern  for  the  "B"  case  causes  the  stage  to  operate  at  a 
better  condition  than  was  experienced  with  uniform  flow. 

This  illustrates  that  the  way  a compressor  responds  to 
radial  distortion  must  depend  on  the  spanwise  work 


92 


AEOC-TR -77-80 


distribution.  Therefore,  the  one-dimensional  stage 
characteristic  used  up  to  this  point  will  not  permit  compu- 
tation of  radial  distortion  effects.  As  noted  in  Fig.  38, 
the  distortion  induced  shift  in  radial  work  distribution,  as 
well  as  the  distorted  flow  pattern  itself,  induces  radial 
crossflows . 

The  three  effects  which  must  be  considered  for 
radial  distortion  analysis  are  then  as  follows. 

Radial  variation  of  compressor  entry  flow 
properties . These  are  automatically  treated  by  the  model 
because  of  the  division  of  the  compressor  into  control 
volumes  radially  as  well  as  circumferentially. 

Radial  variation  of  work  and  force  distribution. 
Stage  characteristics  may  be  specified  for  each  radial 
location  in  the  same  manner  as  a single,  average  stage 
characteristic  was  specified  for  the  one-dimensional 
analysis.  For  the  example  of  Fig.  34,  each  stage  would  have 
a hub,  a midspan,  and  a tip  set  of  stage  characteristics. 

For  the  XC-1  compressor  used  in  example  calcu- 
lations in  Chapter  V,  the  compressor  was  divided  radially 
into  five  annuli  (jn  = 6) . Stage  characteristics  for  each 
annulus  were  not  available,  but  radial  distributions  of 
total  pressure  and  total  temperature  at  each  stage  discharge 
were  given  in  [43]  for  one  operating  point  at  the  corrected 
speed  setting  analyzed.  These  profile  data  were  used  to 
compute  a correction  to  stage  pressure  ratio  and  stage 
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‘tempera'ture  ratio  as  a function  of  radius  for  each  stage; 
thus,  in  effect,  producing  radially  variant  stage  charac- 
teristics. A total  temperature  ratio  correction  factor. 


CF^(r) 


TR(r) 

TR 

avg 


t 


and  a total  pressure  ratio  correction  factor. 


(137) 


CF^(r) 


PR(r) 

PR  ' 
avg 


(138) 


were  determined  for  each  stage  from  the  profile  data  of 
[43] . The  stage  characteristics  as  functions  of  radius  then 


become 


and 


(139) 


Inherent  in  this  correction  is  the  assumption  that  the 
pressure  ratio  and  temperature  ratio  profiles,  normalized  by 
their  averages,  do  not  vary  greatly  over  the  remge  of  flow 
coefficient  of  interest.  Examination  of  data  from  other 
compressors  where  profile  data  were  available  over  a wider 
range  of  operating  conditions  than  used  in  this  study  showed 
the  profiles'  shapes  to  be  consistent  (within  ±5%,  Euid 
usually  better)  up  to  the  point  of  stall. 
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Radial  crossflows.  Racial  crossflows  are  computed 
as  described  previously.  A stage  coefficient  correction  for 
the  distortion-induced  radial  flow  velocity  across  the 
control  volume  face  is  made  by  computing  tj;  as 


Figure  39a  shows  a compressor  subjected  to  a simple  circum- 
ferential distortion  pattern.  Half  of  the  face  of  the 
compressor  is  in  a lower  than  average  total  pressure  region 
and  half  is  above.  The  example  distortion  waveform  is  a 
simple  sine  wave.  Upstream  of  the  compressor  face  a length 
of  one  compressor  diameter  or  more,  the  static  pressure 
variation  across  the  entry  duct  caused  by  the  presence  of 
the  compressor  in  the  distorted  flow  field  will  be  nearly 
zero.  The  total  pressure  distortion  at  that  location  is 
then  principally  an  axial  velocity  distortion.  Depending  on 
design  characteristics  and  operating  conditions  of  the 
compressor,  a portion  of  the  high  velocity  region  will  be 
diffused,  creating  a circumferential  static  pressure 
gradient  at  the  compressor  face.  This  results  in  the 
generation  of  a swirl  velocity  (Fig.  39a)  carrying  mass  from 
the  high  to  the  low  pressure  regions  of  the  compressor  face. 
(This  phenomenon  is  described  in  more  detail  in  [27,  38, 
54].)  Figure  39a  shows  the  influence  of  the  variation  in 
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axial  and  swirl  velocities  at  four  circumferential  locations. 
At  ft  ® 0,  there  is  no  swirl  velocity  and  axial  velocity  is 
high;  therefore,  a is  low  and  the  stage  is  operating  at  a 
lightly  loaded  condition.  At  6 = 90®,  the  axial  velocity  is 
still  fairly  high,  and  the  swirl  component  exists,  and  is  in 
the  direction  of  rotor  rotation.  The  swirl  component 
combines  with  wheel  speed  to  reduce  the  angle  of  attack  at 
this  point.  At  6 = 180®,  the  cixial  velocity  is  at  a 
minimum,  thus  causing  a high  angle  of  attack  and  high,  even 
stalled,  stage  loading.  At  0 = 270®,  the  axial  velocity  is 
increased  once  more,  but  the  swirl  velocity  component  is 
counter  to  rotor  speed,  thereby  causing  a higher  angle  of 
attack  and  stage  loading. 

To  a lesser  extent,  the  same  flow  situation  exists 
at  each  stage  entry.  The  effects  on  the  stage  character- 
istic operating  points  are  accounted  for  by  computation  of 
the  flow  coefficient  as 


= 


u 


U , + V 
wh  sw 


(142) 


The  swirl  velocity,  is  simply  the  circumferential 

velocity  component,  w,  at  that  particular  location. 

Unsteady  effects.  Figure  39b  illustrates  unsteady 
airfoil  effects  for  the  same  circumferential  distortion 
pattern.  As  indicated  in  Fig.  39b,  the  axial  and  swirl 
velocities  caused  the  flow  coefficient  to  decrease,  even  to 
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below  the  steady-state  stage  stall  value.  The  distortion 
pattern  is  stationary.  The  rotor  blades  pass  through  the 
spatially  variant  stationary  pattern  with  each  rotation, 
and,  from  their  point  of  reference,  the  flow  is  cyclically 
unsteady.  For  a sinusoidal  pattern,  the  frequency  felt  by  a 
rotor  blade  is 


LJ  _ N 
2it  ” 60  ' 


(143) 


where  N is  the  mechanical  rotor  speed  in  rpm. 

If  the  distortion  pattern  contained  two  cycles  of 
a sine  wave  instead  of  one  (i.e.,  two  low  pressure  regions), 
the  frequency  would  double.  The  reduced  frequency  for  the 
blade  is 


k 


2irfC 

^rel 


NC 
60  V 


rel 


(144) 


As  discussed  in  Chapter  II,  the  blade  forces  cannot  instan- 
taneously respond  to  rapid  changes  in  angle  of  attack. 
Therefore,  as  a blade  passes  through  the  low  pressure  region 
(region  of  low  * , Fig-  39b),  the  blade  forces,  thus  do 
not  follow  the  <|)  variation  in  a steady-state  manner.  Near 
the  stall  point,  if  the  variation  is  sufficiently  rapid,  the 
lag  is  such  that  may  proceed  into  the  stalled  region,  and 
quickly  retreat,  without  blade  stall  actually  occurring 
135] . The  effect  is  more  pronounced  as  reduced  frequency 
increases . 
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Tlic  unsteady  cascade  aerodynamics  arising  from 
circumferential  distortion  arc  approximated  in  the  model 
following  the  treatment  described  in  Chapter  II,  with  an 
empirical  correction  added  for  the  stalled  blade  region  of 
operation.  In  Figs.  11  and  14,  the  unsteady  lift  magnitude 
is  given  as  a function  of  reduced  frequency.  As  was  shown 
in  Chapter  II,  £q.  (58)  , the  dynamic  response  of  \l)  is  the 
same  as  the  dynamic  response  of  C^.  Therefore,  the  ratio  of 
maximum  dynamic  il»  to  the  quasi-steady-state  value  is  the 
same  as  the  corresponding  lift  coefficient  ratio.  Or, 


U! 


max 

ss 


max 

ss 


DLR  . 


(145) 


(DLR  means  dynamic  lift  ratio.) 

Referring  to  the  upper  left-hand  corner  of 

Fig.  39b,  as  a rotor  blade  begins  at  location  k = 1,  the 

value  of  is  As  it  progresses  to  the  low  pressure 

area,  k = 4,  . The  value  of  calculated  on  a 

steady  flow  basis  in  region  k = 4 is  , but,  is  incorrect 

^ss 

because  of  dynamic  effects.  Let  the  actual  value  be  ij;. 

^dyn 


A correction  may  be  constructed  as  follows. 


gyn 


- V. 


- 


ss 


max 

ss 


DLR  . 


(146) 
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Therefore , 


0.  = DLR  + (1  - DLR)  . (147) 

^dyn  ^ss 

Thus,  the  actual  value  depends  not  only  on  the  steady- 
state  value  computed  at  the  k = 4 location,  but  also  on  the 
blade's  history  as  it  passed  from  location  k = 1. 

This  dynamic  lag  algorithm  can  be  more  generally 
written  so  as  to  be  applied  successively  from  one  sector  to 
the  proceeding  sector  in  the  direction  of  rotor  rotation. 

Or, 

''ijk ' ''ijta  • 

Thus,  a lag  is  introduced  into  the  stage  characteristics. 
The  value  of  DLR  used  was  computed  from 

DLR  = ^ ^ ^ , (149) 

V 1 + (Ek)^ 

as  shown  in  Fig.  14.  The  reduced  frequency,  k,  depended  on 
the  shape  of  the  circumferential  distortion  pattern.  The 
value  of  £ used  followed  the  Goethert-Reddy  analysis,  but 
could  readily  be  varied  based  on  measured  response  of  a 
compressor  to  circumferential  distortion  of  varying  pattern 
shape  to  account  for  additional  viscous  effects.  This 
part  of  the  unsteady  correction  predominates  in  the  blade 
unstalled  region  of  operation. 
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Goethert  and  Reddy  point  out  in  [32]  that  blade 
dynamic  response  should  be  expected  to  be  more  sluggish 
because  of  boundary  layer  effects  not  included  in  their 
inviscid  analysis.  This  will  be  particularly  true  when  the 
blade  is  operating  in  its  stalled  region.  As  already 
pointed  out,  the  blade  may  transiently  operate  at  angles  of 
attack  greater  than  the  steady-state  stall  angle  of  attack 
without  actually  stalling.  Experimental  results  on  the 
degree  to  which  this  phenomenon  may  be  present  for  isolated 
airfoils  were  compiled  by  Melick  [36]  and  are  indicated  in 
Fig.  39c.  The  dynamic  stall  lift  coefficient  increment  is 
correlated  to  the  nondimensional  time  derivative  of  angle  of 
attack  for  several  different  types  of  isolated  airfoils. 

This  effect  was  incorporated  in  the  model  by 
computing  the  derivative  of  angle  of  attack  from  the  change 
in  flow  coefficient  as  a blade  passed  through  the  distortion 
pattern  using  Eg.  (34)  to  relate  flow  coefficient  to  angle 
of  attack.  A stall  lift  coefficient  increment  was  then 
calculated  using  the  empirical  relationship  of  Fig.  39c. 

The  resulting  lift  coefficient  increment  was  then  used  to 
compute  an  equivalent  increment  in  stage  pressure  and 
temperature  coefficients  using  Eq.  (39) , neglecting  the 
influence  of  the  blade  drag  coefficient.  With  this  infor- 
mation, the  stage  characteristic  curve  was  translated 
upward,  at  the  unstalled  stage  characteristic  slope  existing 
near  the  stage  characteristic  peak,  in  a manner  analogous  to 
the  - a curve  translation  of  Fig.  39c.  This  translation. 
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when  combined  with  the  dynamic  lag  effects  of  Eq.  (149) 
yields  the  dynamic  stage  characteristic  behavior  in  the 
stalled  stage  region  depicted  in  Fig.  39b. 

The  Goethert-Reddy  analysis  [32]  and  Fig.  11 
indicated  less  dynamic  effects  with  a cascade  of  airfoils 
than  are  experienced  with  a single,  isolated  airfoil,  it 
might  well  be  anticipated  that  the  stall  region  behavior  of 
the  cascade  would  be  less  than  the  isolated  airfoil  case  as 
well.  For  that  reason,  the  "stall  overshoot"  constant,  m. 
Fig.  39c,  was  chosen  at  four  different  values,  and  influence 
on  the  solution  determined  parametrically.  A value  of  m 
between  one  and  ten  was  found  to  be  reasonable.  Details  of 
the  parametric  investigation  are  discussed  in  Chapter  V. 

IV.  METHOD  OF  SOLUTION 

In  this  section,  the  approximation  to  the  time 
derivative  of  the  volume  integrals  in  the  governing 
equations  is  described.  The  resulting  model  equations, 
including  effects  of  the  simplifications  already  described, 
are  written  and  the  solution  procedure  described. 

Volume  Integral  Time-Derivative  Approximation 

In  each  of  the  governing  equations,  Eqs.  (119), 
(120) , and  (125)  , a partial  time  derivative  of  a volume 
integral  exists . The  general  form  is 
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y ( z,r,s,t  ) d(vol) 


{ z,r,s,t  ) d(vol)  . 


(150) 


It  is  desired  to  approximate  Y in  terms  of  a control  volume 
surface  value  of  y.  By  the  mean  value  theorem,  there  exists 
some  point  (z,r,i)  interior  to  the  control  volume,  such  that 


{ 2,r,s,t  ) d(vol)  = (vol)  |^(z,r,s,t)  . (151) 


Let  be  a value  of  the  integrand  on  either  the  forward  or 
aft  control  volume  face  normal  to  the  z-axis.  The  time 
derivative  at  the  mean  value  point  is  then  related  to  the 
surface  time  derivative  by  expanding  around  the  surface 


point. 


^ [it  [ I'' 

^ [it  [ 

A [It  [ I'' 


(152) 


where 


o < a,  b,  c < 1 . 
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Also, 


and 


The  approximation, 


As  = 


- s 


Y « (vol)  ^ , 


(153) 


(154) 


is  chosen  for  use  in  the  model.  Multiplying  Eg.  (152)  by 
the  volume,  and  rearranging  the  order  of  differentiation 
yields  an  estimate  for  the  truncation  error. 


£ = (vol)  ^ 4s  j . (155) 

The  spatial  derivatives  are  finite;  therefore,  as 

Az  , Ar  , As  0 , 


the  error,  e,  vanishes,  thus  the  approximation  is  consistent 
of  order  one. 

Numerical  stability  was  investigated  by  running 
with  different  time  step  sizes  as  in  the  one-dimensional 
analysis.  A time  step  size  of  approximately  one-half  the 
value  indicated  by  the  Courant-Friedrichs-Lewy  criterion 
was  used . 
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The  actual  approximations  made  were. 
Mass 


3 

Tt 


L 


p d(vol) 


vol . . , 
i]k 


5t  '>ipjk 


Axial  Momentum 


(156) 


/ 


^ ( ( p U ) d (vol)  = vol  ^ ( p U ) , ^^  , (157) 


3t  ' ^ ' ijk 


vol . 
i-:k 


Energy 
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vol.  .. 
i]k 


= vol 


r r . + v^ + 1 1 
[ “"L  ^ j Jipjk 


(158) 


Model  Equations 

Using  the  approximations  and  simplifications  pre- 
viously described  yields  the  equations  comprising  the 
distortion  model. 


Governing  Equations 

Mass,  axial  momentum,  and  energy  were  used.  The 
assumption  of  small  v and  w permits  ignoring  v and  w terms 
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in  some  instances.  Also,  radial  and  circumferential  flows 
were  determined  by  the  crossflow  approximations,  thus 
replacing  the  radial  and  circumferential  momentum  equations. 
The  resulting  equations  are  as  follows. 

Mass 


8t 


- WZ 


ipjk 


Axial  Momentum 


[WZ.  ..  + WR.  ..  + WC.  .. 
i]k  ijk  ijk 

- WR.  . . - WC.  .,  1 . 

rjpk  i]kp  J 


(159) 


a ( P U ) 


ijk  _ 


vol 


ijk 


F.  ..  + IMP.  ..  - IMP 
i]k  i]k 


■ 

ipjk^ 


(160) 


In  Eq.  (160)  the  contribution  of  the  pressure  integrated 
over  the  projected  areas  of  the  radially  oriented  surfaces 
is  combined  with  the  force  term,  FZ,  of  Eq.  (120)  to  pro- 
duce F.  The  force,  F,  is  obtained  from  stage  characteristic 
information,  as  was  done  in  the  one-dimensional  case.  The 
impulse  terms  are 

IMP  = PS  AZ  + WZ  U . (161) 


The  momentum  convected  across  the  radial  and  circumferential 
boundaries  is  negligible  compared  to  the  axial  impulse  terms. 
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Energy 


ax. 


1 

(vol) 


ws 


ijk 


+ Q441,  + HZ 


ijk  '^ijk 


ijk 


-“ipllc 


(162) 


In  Eq.  (162)  , 


X = p 


(163) 


2 2 

The  V and  w terms  in  £q.  (125)  are  negligible  compared  to 

U2. 


The  terms  neglected  In  Eqs.  (160)  and  (162)  are 
second-order  terms  determined  by  nondimensionalization  and 
order  of  magnitude  analysis  of  the  equations.  The  terms 
were  ignored  because  their  influence  on  the  solution  was 
small.  Also,  inclusion  of  terms  caused  the  numerical  solu- 
tions to  be  generally  less  stable  and  to  oscillate  with 
small  amplitude  when  a steady-state  condition  was  reached. 


Solution  Procedure 

The  solution  procedure  used  on  the  digital 
computer  paralleled  the  approach  used  for  the  one- 
dimensional model  outlined  in  Fig.  16.  Initial  values  of 
the  dependent  variables,  p,  p U,  and  X,  were  specified  for 
every  control  volume  location,  ijk,  to  start  the  solution. 
For  simplicity,  the  solution  was  always  started  from  a 
steady,  uniform  flow  condition,  thus  easing  the  chore  of 
computing  initial  conditions. 
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Forward  boundary  conditions.  Forward  boundary 
conditions  were  total  pressure  and  total  temperature, 

, and  T-  . , , specified  at  each  control  volume  face. 

± ,J  f K.  1 ,2  fK 

Entry  plane  conditions  were  then  computed  using  the  same 
relationships  as  indicated  in  Fig.  16,  but  for  every  control 
volume  face  instead  of  just  one.  Thus,  the  subscript,  1,  in 
Fig.  16,  entry  plane  conditions,  is  replaced  with  subscript, 
1, j ,k. 


Interior  plane  values.  The  same  relationships 
indicated  in  Fig.  16  are  used  again.  Once  more,  the 
subscript,  i,  must  be  replaced  with  i,j,k,  and  the  quanti- 
ties computed  at  each  control  volume  face. 

Crossflow  values.  Crossflows  among  the  control 
volumes  are  computed  using  Eqs.  (127)  through  (136)  for  each 
control  volume.  For  all  calculations  made  in  this  study, 
crossflow  across  the  compressor  walls  was  zero. 

Exit  plane  conditions.  Only  unchoked  flow  at  the 
exit  was  used  for  this  study;  thus,  a uniform  static 
pressure  was  specified  at  every  control  volume  aft  face  at 
the  compressor  exit  plane.  Exit  plane  calculations  follow- 
ing the  unchoked  portion  of  the  one-dimensional  procedure 
(Fig.  16,  exit  plane  conditions)  were  therefore  used  with 
the  index,  in,  replaced  by  in,j,k,  and  computations  made  for 
all  control  volumes. 
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Force  and  work  computations.  The  same  procedure 
used  in  the  one-dimensional  procedure  was  used  again,  but 
individually  at  each  control  volume.  Local  values  of  flow 
coefficient  were  computed  and  adjusted  for  distortion 
effects  according  to  Eqs.  (141)  and  (142) . Further  cor- 
rections to  stage  characteristics  for  radial  variations  and 
for  unsteady  effects  were  also  made. 

Time-derivative  computations.  The  derivatives  of 
the  dependent  variables,  o,  p U,  and  X,  are  computed  from 
Eqs.  (159),  (160),  and  (162)  for  every  control  volume.  The 
values  are  then  integrated  to  provide  values  of  the  vari- 
ables at  the  next  time  step.  The  sequence  is  then  repeated, 
with  the  boundary  conditions  changing  in  accordance  with  the 
specified  problem  under  investigation  and  the  computations 
continued  until  the  desired  solution  is  obtained. 

Additional  numerical  stability  considerations.  To 
prevent  numerical  instabilities  in  the  model  solution  caused 
by  error  buildup  in  the  crossflow  calculations  when  the 
driving  pressure  differences  from  control  volume  to  control 
volume  are  very  small,  the  pressure  differences  are  tested 
to  assure  a significant  pressure  gradient  exists.  If  a 
significant  difference  exists,  crossflow  is  calculated;  if 
not,  crossflow  is  set  to  zero. 
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CHAPTER  V 

EXAMPLE  PROBLEMS  WITH  COMPARISON  TO  EXPERIMENTS, 
DISTORTION  MODEL 

One  steady-state  combined  radial  and  circumfer- 
ential distortion  case  for  the  XC-1  compressor  is  given. 
Detailed  experimental  data  were  available  for  this  case  and 
comparisons  between  model  computations  and  experimental 
measurements  are  presented.  An  example  of  time-variant 
distortion  for  the  XC-1  compressor  is  given.  Pure  radial 
and  pure  circumferential  distortion  effects  are  computed. 

A pure  circumferential  temperature  distortion  case  is  also 
presented.  Specific  experimental  data  were  not  available 
for  the  time-variant  distortion,  pure  radial  and  pure 
circumferential  distortion  patterns  for  the  XC-1  compressor. 
Data  from  other  compressor  tests  were  used,  therefore,  to 
check  the  model  results  for  proper  trend  and  order  of 
magnitude . 

In  all  cases  discussed  in  this  chapter,  the 
nominal  compressor  entry  conditions  were  one  atmosphere 
pressure  and  518. 7®R  temperature. 

I.  COMBINED  DISTORTION  PATTERN,  XC-1  COMPRESSOR 

Figure  40  shows  a contour  map  of  the  combined 
circumferential  and  radial  total  pressure  distortion  pattern 
to  which  the  XC-1  compressor  was  subjected.  The  pattern  was 
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generated  using  a distortion  screen  to  simulate  a distortion 
pattern  present  in  the  inlet  of  a VTOL  aircraft  [43] . The 
pattern  of  Pig.  40  is  for  a corrected  rotor  speed  of  90%. 

The  pattern  amplitude  was  10.2%.  The  lowest  pressure  region 
was  centered  at  approximately  210  degrees,  in  the  tip  region. 

The  distortion  pattern  was  input  to  the  model  by 
dividing  the  compressor  annulus  as  indicated  in  Fig.  41-. 

Each  stage  was  divided  into  30  control  volumes.  The  total 
pressure  at  the  compressor  entry  plane  was  averaged  over 
each  control  volume  face  to  obtain  the  values  indicated  in 
Fig.  41.  The  averaging  resulted  in  reduction  in  the 
apparent  distortion  pattern  amplitude  of  approximately  1.0%. 
The  compressor  was  divided  axially  into  four  stages  (five 
axial  locations)  as  shown  in  Fig.  34.  Inlet  total  tempera- 
ture was  constant. 

The  compressor  discharge  static  pressure  was 
specified  as  the  aft  boundary  condition  and  was  ramped 
upward  to  load  the  compressor  to  the  stability  limit.  The 
process  was  conducted  with  and  without  distortion.  The 
results  are  shown  in  Fig.  42,  and  are  compared  to  experi- 
mental measurements.  A reduction  of  5.3%  in  stability  limit 
pressure  ratio  caused  by  distortion  was  computed.  The 
results  agreed  well  with  experimental  measurements. 

Figure  43a  presents  a detailed  comparison  of  the 
circumferential  total  pressure  profiles  in  the  tip  region  of 
the  compressor  at  each  stage  entry  and  exit.  The  distortion 
pattern  is  seen  to  attenuate  until , at  the  compressor 
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discharge,  the  variation  is  within  1%  of  the  average  value. 
The  agreement  between  the  model  and  experiment  is  reasonably 
good,  usually  within  1%,  but  occasionally  different  by  as 
much  as  2%. 

Two  mechanisms  serve  to  attenuate  the  circumfer- 
ential distortion.  The  first  stems  from  the  fact  that  in 
the  low  pressure  region,  the  flow  coefficient  is  lower; 
thus,  if  the  stage  is  not  stalled,  it  will  work  harder 
{produce  a higher  pressure  and  temperature  coefficient)  than 
the  high  pressure  regions.  This  effect  tends  to  drive  the 
stage  discharge  pressures  of  the  low  and  high  pressure 
regions  toward  equality.  Figure  43b  shows  the  corresponding 
total  temperature  profiles  and  indicates  the  transposition 
of  the  pressure  distortion  to  temperature  distortion.  The 
second  mechanism  working  to  smooth  out  the  circumferential 
pressure  profile  is  circumferential  crossflow,  which  was 
discussed  in  Chapter  IV.  The  slight  phase  shift  between  the 
experimental  and  computed  profiles  may  be  caused  in  part  by 
slightly  different  crossflow  in  the  actual  compressor  than 
that  allowed  for  in  the  model.  (Crossflow  effects  are  dis- 
cussed in  more  detail  later  in  this  chapter.) 

Another  reason  for  the  apparent  phase  shift  is  the 
finite  size  of  the  control  volumes.  A finer  circumferential 
division  would  reduce  any  skew  of  the  computed  profile 
waveform  caused  by  averaging.  The  profile  waveform  ampli- 
tude is  also  reduced  by  the  averaging  effect  of  the  finite 
control  volume,  just  as  occurred  between  Figs.  40  and  41. 
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Thus,  the  model  profile  amplitudes  should  be  expected  to  be 
1%  or  so  less  than  the  experimentally  observed  profiles. 

Figure  44  shows  the  attenuation  of  the  distortion 
pattern  through  the  compressor  for  the  tip  (Fig.  44a] , 
mid-span  (Fig.  44b),  and  hub  (Fig.  44c)  regions.  Agreement 
between  the  computed  and  experimental  results  is  generally 
within  l?j  when  allowance  is  made  for  the  averaging  effects 
of  the  finite  control  volumes  of  the  model.  Pressure  dis- 
tortion reduction  and  temperature  distortion  production  are 
similar  for  the  tip  and  mid-span  regions.  In  the  hub 
region,  however,  both  the  computed  and  experimental  measure- 
ments indicate  amplification  of  the  distortion  across  the 
first  stage.  This  occurs  because  of  the  shallower  stage 
characteristics  which  exist  in  the  hub  region.  The  agree- 
ment is  not  as  good  in  the  hub  region  as  in  the  tip  and 
mid-span  regions,  indicating  some  inaccuracy  in  the 
estimation  of  the  hub  region  stage  characteristics. 

Figure  45  presents  a comparison  of  radial  total 
pressure  and  temperature  profiles  in  the  region  of  maximum 
radial  distortion  (at  210  degrees) . Figure  45a  shows  agree- 
ment within  1 to  2%  except  in  the  hub  region  at  the  second 
stage  entry.  Once  more,  this  indicates  some  inaccuracy  in 
the  hub  region  stage  characteristic.  Figure  45b  shows  the 
total  temperature  radial  profiles.  Agreement  within  1%  was 
generally  achieved. 

Figure  46  presents  similar  profiles  for  the 
compressor  operating  without  distortion.  Agreement  between 
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computed  and  experimental  results  was  generally  within  1%. 
Comparing  Figs.  45  and  46  indicates  that  the  radial  dis- 
tortion with  low  pressure  in  the  hub  decays  rapidly  in  the 
compressor,  but  some  modification  to  the  profiles  does 
occur. 

Computation  of  proper  distributions  within  the 
compressor  is  important  to  stability  limit  determination 
because  it  determines  the  local  stage  flow  coefficient 
distribution  which  in  turn  controls  stability.  Figure  47 
presents  the  approximate  flow  coefficient  distribution 
computed  for  each  stage  just  prior  to  reaching  the  com- 
pressor stability  limit.  The  flow  coefficient  values  are 
normalized  by  the  stage  stall  value.  As  should  be  expected, 
the  distribution  approximates  the  shape  of  the  total 
pressure  distortion  pattern.  Stage  1,  Fig.  47a,  is  almost 
stalled  in  a small  region  at  the  tip.  Stage  2,  Fig.  47b,  is 
actually  operating  with  the  tip  stalled  in  the  distorted 
region.  Figure  47c  shows  a similar  pattern  for  stage  3,  but 
it  is  not  stalled.  Stage  4,  Fig.  47d,  is  stalled  over  a 
wide  arc  from  hub  to  tip,  and,  along  with  stage  2,  is  indi- 
cated by  the  model  to  cause  the  flow  breakdown  leading  to 
total  compressor  instability. 

The  stalled  region  in  stage  4 is  seen  to  extend 
beyond  the  concentrated  low  pressure  region  of  the  imposed 
distortion  pattern.  The  "pattern  spreading"  is  caused  by 
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the  attenuation  mechanisms  discussed  previously.  When  a 
sufficiently  large  region  becomes  affected,  total  insta- 
bility will  occur. 

Control  Volume  Size  Selection 

The  example  problem  discussed  in  this  section 
illustrates  the  division  of  the  compressor  into  finite 
control  volumes.  The  circumferential  division  (60-degree 
arc  sectors)  was  made  because,  as  Goethert  and  Reddy  point 
out  in  [32] , it  is  at  approximately  60  degrees  that  each  arc 
sector  acts  very  much  as  a separate  compressor  (as  parallel 
compressor  theory  would  suggest) ; thus,  the  effects  of  the 
unsteady  cascade  aerodynamic  and  crossflow  approximations 
are  less  crucial  to  achieving  an  accurate  representation  of 
the  compressor  with  regard  to  circumferential  distortion. 

As  already  noted,  some  error  is  introduced  into  the  model 
results  because  of  local  distortion  pattern  variations 
within  a 60-degree  arc  sector  which  are  averaged  out  by  the 
approximating  process. 

The  radial  division  selection  (five  rings)  was 
made  based  on  the  number  of  radial  locations  at  which 
experimental  total  temperature  and  total  pressure  profile 
data  were  available  for  approximation  of  the  radial  vari- 
ation of  stage  characteristics.  In  general,  the  selection 
of  the  radial  division  must  depend  on  the  strength  of  the 
radial  gradients  of  the  distortion  pattern  and  the  internal 
radial  profiles  within  the  compressor.  Divisions  which 
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hold  avoraqing  error  to  within  1 or  2%  should  bo  quite 
sufficient.  Also,  the  radial  crossflow  representation  in 
the  model  is  relatively  crude.  Therefore,  it  should  not  be 
expected  that  finer  and  finer  radial  division  of  the  com- 
pressor would  yield  greater  and  greater  accuracy  of  the 
computed  results.  Additionally,  the  stage  characteristic 
data  are  based  on  one-dimensional  (or,  at  best,  averaged 
over  a radial  distance)  data,  thus,  a very  fine  division 
radially  would  not  be  consistent  with  the  resolution  of  the 
input  empirical  information. 

II.  TIME-VARIANT  DISTORTION 

Figure  48  illustrates  time-variant  distortion. 

The  total  pressure  at  the  engine  entry  plane  fluctuates  with 
time  and  also  varies  spatially  across  the  compressor  face. 
The  fluctuations  are  caused  by  unsteady  flow  processes  in 
the  inlet,  such  as  unsteady  shock  wave-boundary  layer 
separation  [6,  13,  15).  The  unsteady,  nonuniform  flow 
causes  distortion  patterns  at  the  engine  entry  which  vary 
with  time,  as  illustrated  in  Fig.  48.  The  spatially 
averaged  pressure  also  oscillates. 

The  XC-1  compressor  model  was  subjected  to  time- 
variant  distortion  by  superimposing  a time-dependent 
waveform  on  the  distortion  pattern  of  Fig.  40.  Figure  49 
shows  the  frequency  spectrum  of  the  superimposed  fluctu- 
ations. The  spectrum  is  representative  of  those  measured 
experimentally  in  aircraft  inlets  [15].  The  waveform  was 
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applied  in  two  parts.  One  part  was  phase  shifted  over  the 
compressor  face  so  as  to  alternately  stretch  and  shrink  the 
steady-state  distortion  pattern.  The  other  part  was  added 
uniformly  to  produce  timewise  fluctuation  in  the  spatial 
average.  (Details  of  time-variant  distortion  pattern 
analysis  and  synthesis  are  discussed  in  [42].) 

Two  cases  were  computed  for  the  XC-1  compressor 
with  time-variant  distortion.  One  was  for  an  overall  root- 
mean-square  (rms)  amplitude  of  2%  and  the  other  was  for  4%. 
With  these  inlet  flow  fields  present,  the  compressor  exit 
pressure  was  ramped  upward  until  instability  occured. 

Figure  50a  shows  the  loss  in  stability  limit  pressure  ratio 
correlated  with  the  rms  amplitude.  A reduction  results,  as 
should  be  expected,  and  is  of  similar  magnitude  as  the 
experimental  results  of  [7]  for  a different  compressor. 

The  model  computed  results  were  also  correlated  to 
the  aunplitude  of  the  instantaneous  distortion  pattern 
existing  just  prior  to  compressor  instability,  as  shown  in 
Fig.  50b.  All  pattern  shapes  were  similar,  thereby  per- 
mitting use  of  the  simple  distortion  pattern  magnitude 
descriptor.  The  results  generalize  well  with  steady-state 
distortion  results,  thus  providing  a degree  of  analytical 
conformation  to  the  quasi-steady-state  distortion  approach 
currently  used  in  experimental  stability  analysis.  (See 
Chapter  I and  [7]  for  further  details  on  that  approach.) 

The  results  shown  in  Fig.  50  should  not  be  sur- 
prising because  the  XC-1  compressor  is  very  small  (about 
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17  inches  in  diameter  and  8 inches  long)  and  has  high  rotor 
speed  (about  16,230  rpm  at  90%  speed).  Therefore,  fluctu- 
ations in  the  frequency  range  considered  are  sufficiently 
low  to  permit  the  instantaneous  distortion  pattern  to  dwell 
at  the  entry,  and  penetrate  through  the  compressor  in  a 
quasi-steady  manner.  For  example,  in  this  case,  approxi- 
mately 0.0037  second  is  required  for  one  rotor  revolution, 
and  0.0017  second  is  required  to  convect  a particle  through 
the  compressor  (average  axial  velocity  of  400  ft/sec) . The 
period  of  the  highest  frequency  component  in  the  spectrum  is 
0.0038  second.  Thus,  the  distortion  effects  will  be  much 
like  steady-state  distortion. 

If  larger  cunplitude  timewise  fluctuations  are  con- 
sidered, particularly  if  the  compressor  is  combined  with 
upstream  ducting  and  a combustor,  planar  oscillation  effects 
similar  to  those  discussed  in  Chapter  III  would  show  an 
effect  also. 


III.  PURE  RADIAL  AND  CIRCUMFERENTIAL 
DISTORTION  PATTERNS 

Three  additional  types  of  distortion  patterns  were 
computed  to  gain  some  insight  into  the  influence  of  pattern 
g0yg]fity  and  shape  variations.  The  patterns  considered  arc 
pure  radial  pressure  distortion,  pure  circumferential 
pressure  distortion,  and  pure  circumferential  temperature  dis- 
tortion. No  specific  experimental  data  for  the  XC-1 
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compressor  were  available  for  these  cases.  Comparisons  to 
experimental  results  from  other  compressors  are  made 
therefore. 

Radial  Pressure  Distortion 

Figure  51a  illustrates  the  results  of  loading  the 
XC-1  compressor  model  to  the  stability  limit  with  tip  radial 
distortion  (low  pressure  in  the  tip  region) . Pattern 
severities  of  5,  10,  and  20%  were  considered.  Because  of 
the  relatively  high  tip  loading  of  this  compressor,  Fig.  46, 
it  appears  to  be  more  sensitive  to  tip  radial  distortion. 
Figure  51b  shows  the  results  for  the  hub  radial  case  and 
indicates  the  compressor  to  be  more  tolerant.  Figure  52 
shows  comparisons  of  the  XC-1  model  results  with  experi- 
mental data  from  a J85-13  test  [16] . Similar  results  were 
obtained,  except  that  the  J85-13  was  even  more  tolerant  of 
hub  radial  distortion  than  the  XC-1  was  computed  to  be. 

Circumfere.ntial  Pressure  Distortion 

Figure  53  shows  the  effect  of  increasing  the 
severity  of  a 60-degree-arc  circumferential  distortion 
pattern.  Pattern  magnitudes  of  5,  10,  15,  and  20%  were 
computed.  Increasing  loss  of  stability  limit  pressure  ratio 
and  airflow  was  computed.  Figure  54  shows  the  model 
results  compared  to  a range  of  experimental  data  from  other 
compressors.  The  model  results  appear  reasonable. 

Figure  55a  depicts  the  variation  of  airflow  in  the 
distorted  region.  Crossflow  into  the  distorted  region 
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rapidly  increases  the  flow  in  that  sector.  Also  shown  in 
Fig.  55a  is  the  crossflow  into  an  approximately  90-degree 
distorted  region  computed  for  the  XC-1  compressor  from 
experimental  measurements  [43] . 

Two  different  crossflow  distribution  cases  were 
computed  to  check  the  influence  of  crossflow  on  computed 
compressor  behavior.  In  one  case,  the  first  stage  cross- 
flow  was  increased.  The  increase  was  accomplished  by 
increasing  the  first  stage  crossflow  area.  Figure  55b  indi- 
cates the  effect  of  the  crossflow  change  on  compressor 
stability.  The  increased  first  stage  crossflow  case 
matched  the  experimental  distribution  better,  and  caused  a 
higher  stability  limit  pressure  ratio  with  distortion 
present.  The  improved  performance  comes  about  by  allowing 
crossflow  to  relieve  downstream  stage  loading  in  the  dis- 
torted region.  The  crossflow  distribution  more  closely 
matching  the  experimental  case  was  used  in  the  previously 
discussed  combined  distortion  case. 

Figure  56  shows  the  influence  on  stability  of 
varying  the  angular  extent  of  the  circumferential  low 
pressure  area  while  holding  the  distortion  pattern  magnitude 
constant.  Three  different  cases,  6 = 30,  60,  and  180 

degrees  are  shown.  The  30-degree  case  has  the  least  loss  in 
stability.  The  60-  and  180-degree  cases  have  approximately 
the  same  percentage  loss  in  stability  limit  pressure  ratio 
(taken  at  constant  corrected  airflow  through  each  stability 
limit  point).  However,  the  flow  and  pressure  ratio 
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reduction  for  the  180-degree  case  was  more  than  that  of  the 
60-degree  case.  Figure  37  shows  the  loss  in  stability  limit 
pressure  ratio  as  a function  of  distortion  pattern  angular 
extent.  Computations  were  made  for  four  different  values  of 
dynamic  lift  coefficient  "stall  overshoot"  constant,  m. 

Fig.  39c.  Also  shown  in  Pig.  57  are  experimental  data  from 
a J85-13  compressor  subjected  to  distortion  of  varying 
angular  extent.  The  experimental  results  fall  between  the 
m = 1 and  m = 10  model  computation  results.  The  m = 40 
results,  corresponding  to  the  center  of  the  isolated  airfoil 
correlation  curve,  Fig.  39c,  produce  virtual  insensitivity 
to  circumferential  distortion  up  to  60  degrees,  which  is  an 
unrealistic  result.  A value  of  m between  1 and  10  is  a 
more  realistic  value.  The  results  shown  in  Fig.  56  are  for 
m = 10 . 

The  analysis  presented  in  Fig.  57  indicates  that 
the  unsteady  cascade  aerodynamic  effects  predominate  in 
determining  the  response  of  the  compressor  to  circumfer- 
ential distortion  below  approximately  60  degrees  angular 
extent. 

Circumferential  Temperature  Distortion 

Figure  58  shows  the  effect  of  varying  the  severity 
of  total  temperature  distortion  on  compressor  stability. 
Temperature  distortion  destabilizes  the  compressor  in  much 
the  same  manner  as  does  pressure  distortion.  The  distorted 
region,  which  is  hotter  than  the  rest  of  the  compressor. 
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Operates  at  reduced  airflow,  thus  at  reduced  flow  coef- 
ficients compared  to  the  rest  of  the  compressor.  Therefore, 
as  the  load  on  the  compressor  is  increased,  the  hot  region 
stalls  first,  eventually  causing  the  entire  compressor  to 
become  unstable. 

Mechanical  rotor  speed  was  constant  during  the 
loading.  Corrected  rotor  speed,  indicated  in  Fig.  58, 
decreased  as  much  as  1.1%  because  of  the  increased  average 
temperature  over  the  face  of  the  compressor . 

Little  influence  of  crossflow  was  found,  indi- 
cating that  basic  parallel  compressor  theory  (with 
appropriate  unsteady  cascade  corrections)  should  give 
reasonable  predictions  for  temperature  distortion. 

Figure  59  shows  a comparison  of  model  computed 
stability  limit  pressure  ratio  reduction  to  experimental 
results  from  other  compressors.  The  predicted  trend  is 
reasonable,  although  more  pressure  ratio  reduction  is 
computed  than  was  experienced  on  the  other  compressors. 
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CHAPTER  VI 

SUMMARY  OF  RESULTS  AND  RECOMMENDATIONS 
FOR  FURTHER  WORK 

The  overall  objective  of  this  work  was  to  develop 
an  improved  analytical  method  for  the  computation  of  turbine 
engine  compressor  stability  loss  caused  by  a wide  range  of 
different  types  of  external  flow  disturbances.  Another 
objective  was  to  gain  understanding  of  the  physical 
processes  leading  up  to  instability  as  induced  by  external 
disturbances.  The  understanding  was  provided  by  the 
development  of  an  analytical  model  and  by  application  of  the 
model  to  example  problems  with  comparisons  to  experimental 
results.  The  objectives  were  accomplished  in  two  major 
steps.  In  the  first  step,  a one-dimensional,  time-dependent 
math  model  was  developed  and  applied  to  a variety  of  steady 
and  time-dependent  planar  disturbances  which  are  destabi- 
lizing to  the  compressor.  In  the  second  step,  the  concepts 
of  the  one-dimensional  model  were  extended  to  three 
dimensions  to  allow  the  development  of  a time-dependent 
model  for  computation  of  radial,  circumferential  and 
combined  distortion  effects  on  stability.  The  distortion 
model  was  then  applied  to  the  computation  of  distortion 
cases  with  comparisons  to  experimental  results.  The  results 
of  those  two  steps  are  summarized  as  follows. 
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I.  ONE-DIMENSIONAL,  TIME-DEPENDENT  ANALYSIS 
Model  Description 

The  one-dimensional  model  satisfies  the  time- 
dependent  mass,  momentum  and  energy  equations.  The  full, 
nonlinear  form  of  the  equations,  not  used  in  earlier  works, 
was  used.  The  compressor  is  divided  into  control  volumes  on 
a stage-by-stage  basis.  The  force  and  shaft  work  applied  to 
the  fluid  is  determined  from  empirical  steady-state  stage 
characteristics  modified  for  unsteady  cascade  airfoil 
effects.  The  solution  of  the  model  governing  equations  is 
obtained  through  finite  difference  approximation  of  the 
governing  differential  equations.  The  resulting  system  of 
ec[uations  was  solved  on  a digital  computer. 

Example  Problems 

Steady-state  stability  limit.  Three  different 
compressors  were  modeled.  They  were,  an  Allison  XC-1  lift 
engine  four-stage  compressor,  a NACA  eight-stage  research 
compressor,  and  the  General  Electric  J85-13  turbojet  engine 
compressor.  The  model  was  shown  to  be  able  to  generally 
compute  the  steady-state  stability  limit  pressure  ratio  and 
corrected  airflow  to  within  1 to  2%  of  experimentally 
determined  values.  Compressor  internal  steady-state  flow 
conditions  in  agreement  with  experimental  values  were  also 
shown  to  be  computed  by  the  model . The  initial  flow 
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breakdown  process  in  the  compressor  at  the  stability  limit 
was  shown  to  be  reasonably  well  reproduced  by  the  model . 

Dynamic  stability  limits  (oscillating  inflow 
loading) , NACA  eight-stage  compressor.  The  NACA  eight-stage 
compressor  model  was  driven  to  its  stability  limit  by  planar 
oscillation  of  the  compressor  inlet  pressure  at  frequencies 
from  20  to  320  Hz.  The  model  results  indicated  that,  as 
frequency  increased,  the  amplitude  required  to  cause  insta- 
bility decreased.  It  also  showed  that  the  individual  stages 
transiently  enter  and  retreat  from  their  stalled  regions 
during  the  oscillations  without  total  compressor  instability 
occurring.  As  the  oscillation  amplitude  increased,  however, 
the  time  period  during  which  the  stage  is  in  the  stalled 
region  is  increased,  and  a point  is  eventually  reached  where 
total  compressor  instability  occurs.  The  analysis  showed 
further  that,  for  planar  oscillations  in  the  frequency  range 
analyzed,  the  influence  on  stability  of  the  unsteady  cascade 
correction  to  the  steady-state  stage  characteristics  was 
small . 

The  ability  of  the  model  to  compute  dynamic 
stability  limits  was  shown  to  depend  strongly  on  the  model’s 
ability  to  properly  compute  the  dynamic  response  of  the 
compressor  to  imposed  disturbances . 

Oscillating  inflow  dynamic  response,  J85-13 
compressor . The  unstalled,  dynamic  response  to  planar  entry 
pressure  oscillations  in  terms  of  amplitude  ratio  and  phase 
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angle  were  computed  and  compared  to  experimental  data.  The 
experimental  data  were  obtained  by  NASA  in  a wind  tunnel  at 
Mach  2 . 5 with  the  engine  operating  behind  a supersonic 
inlet.  The  compressor  entry  pressure  oscillations  were  pro- 
duced by  forced  oscillations  of  the  inlet  bypass  doors. 

Data  were  acquired  from  1 to  50  Hz.  Amplitude  ratios  and 
phase  angles  computed  by  the  model  agreed  with  experimental 
values  to  within  10%  and  five  degrees,  respectively.  An 
additional  set  of  experimental  oscillating  entry  pressure 
data  obtained  by  NASA  was  compared  to  the  inlet  data  and 
model  computations.  The  additional  set  was  obtained  in  a 
test  cell  with  an  array  of  small  airjets  blowing  against  the 
primary  engine  flow  in  an  oscillatory  manner  to  produce  the 
pulsations.  The  results  differed  markedly  from  the  inlet- 
engine  results,  probably  because  of  differences  in  the 
boundaries  of  the  dynamic  system  under  consideration. 

Oscillating  compressor  discharge  pressure,  J85-13 
compressor.  The  unstalled,  dynamic  response  of  the  J85-13 
compressor  to  oscillations  in  discharge  pressure  caused  by 
engine  fuel  flow  oscillations  was  computed  using  the  model 
and  compared  to  experimental  results.  A frequency  range  of 
1 to  60  Hz  was  considered.  Agreement  between  model  compu- 
tations and  experimental  results  was  10  to  20%  and  5 to  10 
degrees  on  amplitude  ratio  and  phase  angle,  respectively,  at 
40  Hz.  The  experiment  was  conducted  in  a wind  tunnel  at 
Mach  2.5  with  the  engine  installed  behind  a supersonic  inlet. 
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Three  dynamic  disturbance  paths  were  discovered  to 
exist  at  frequencies  below  20  Hz  and  were  included  in  the 
model  analysis.  The  three  paths  are:  (1)  the  direct  up- 
stream propagation  of  the  pressure  disturbances,  (2)  oscil- 
lation of  stage  pumping  because  of  rotor  speed  oscillations, 
and  (3)  oscillation  of  compressor  entry  pressure  caused  by 
the  supersonic  inlet  shock  wave  oscillations  resulting  from 
the  engine  airflow  oscillation.  The  exeunple  illustrated  the 
influence  of  multiple  disturbance  paths  on  dynamic  behavior 
and  stability. 

Rapid  inlet  temperature  ramp,  NACA  eight-stage 
conpressor . The  NACA  eight-stage  compressor  model  was 
driven  to  instability  by  ramping  inlet  temperature  at  a rate 
of  3600®R  per  second  to  simulate  hot  gas  ingestion  from 
rocket  or  gun  fire.  Ramp  heights  of  100“R  or  over  caused 
instability.  The  results  were  compared  to  a band  of  experi- 
mental data  from  other  compressors.  The  model  computations 
were  within  the  band. 

II.  DISTORTION  ANALYSIS 


Model  Description 

The  concepts  of  the  one-dimensional  model  were 
extended  to  allow  computation  of  the  effects  on  stability  of 
distorted  flow  at  the  compressor  entry  plane.  The  com- 
pressor was  divided  radially  and  circumferentially  as  well 
as  axially  into  control  volumes,  and  crossflow  between  the 
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control  volumes  caused  by  distortion  was  approximated.  The 
governing  equations  were  developed  from  the  three- 
dimensional  , time— dependent  mass,  momentum  and  energy  equa- 
tions for  a finite  control  volume.  Empirical  crossflow 
relationships  were  used  to  replace  the  radial  and  circumfer- 
ential momentum  equations  for  purposes  of  simplification. 

The  radial  work  distribution  variation  effect  necessary  for 
treatment  of  radial  distortion  was  approximated  using 
average  stage  characteristics  modified  as  a function  of 
radius  by  an  empirical  total  pressure  and  total  temperature 
profile  correction.  Crossflow  and  unsteady  cascade  effects 
for  circumferential  distortion  were  built  into  the  model 
also.  The  resulting  set  of  equations  was  solved  using  a 
digital  computer. 

Example  Problems 

Steady-state  combined  radial  and  circumferential 
distortion . The  four-stage,  XC-1  compressor  was  modeled  for 
distortion  by  dividing  each  stage  into  five  radial  and  six 
circumferential  divisions.  This  produced  30  control  volumes 
per  stage,  or  120  control  volumes  for  the  entire  compressor. 
The  model  was  subjected  to  a combined  distortion  pattern  and 
loaded  to  the  stability  limit  using  compressor  discharge 
static  pressure  as  the  aft  boundary  condition.  Detailed 
experimental  data  were  available  for  this  case  and  compari- 
sons were  made.  The  stability  limit  pressure  ratio  and 
corrected  airflow  were  computed  to  within  1,3  to  0.8%  of  the 
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experimental  values,  respectively.  Detailed  circumferential 
and  radial  profile  comparisons  of  total  pressure  and  total 
temperature  were  made.  Agreement  with  experimental  results 
was  reasonable. 

Time-variant  distortion.  A timewise  fluctuating 
total  pressure  was  superimposed  on  the  steady-state  combined 
radial  and  circumferential  distortion  pattern  to  which  the 
XC-1  compressor  model  was  subjected  to  simulate  time-variant 
distortion.  The  reduction  in  stability  limit  pressure  ratio 
was  shown  to  correlate  well  with  the  severity  of  the  instan- 
taneous distortion  pattern  existing  just  prior  to  compressor 
instability.  The  correlation  also  agreed  with  results 
obtained  by  imposing  steady-state  distortion  of  increased 
severity  on  the  model.  This  result  provides  an  analytical 
confirmation  of  the  guasi-steady-state  distortion  method  of 
analysis  currently  used  in  experimental  stability  test  and 
analysis . 

Pure  Radial  and  Circumferential  Distortion 

The  XC-1  compressor  model  was  subjected  to  pure 
radial  pressure  distortion,  pure  circumferential  pressure 
distortion,  and  pure  circumferential  temperature  distortion 
patterns  of  varying  severity  and  their  influences  on 
stability  were  computed.  No  specific  experimental  data  for 
the  XC-1  compressor  subjected  to  these  specific  distortion 
patterns  were  available.  Therefore,  comparisons  were  made 
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to  experimental  results  from  other  compressors.  Trends  and 
general  magnitudes  computed  by  the  model  were  in  agreement 
with  the  experimental  results. 

The  response  of  the  compressor  to  pure  radial 
distortion  was  indicated  by  the  model  analysis  to  depend 
most  strongly  on  the  radial  work  distribution  (radial  stage 
characteristic  variations)  built  into  the  compressor. 
Response  to  pure  circumferential  pressure  distortion  and 
attenuation  of  the  distortion  pattern  through  the  compressor 
depended  on  a combination  of  stage  characteristic  curve 
slope  (steepness) , and  mildly  on  crossflow  distribution. 
Response  to  circumferential  distortion  patterns  with  low 
pressure  region  angular  extent  less  than  60  degrees  was 
indicated  to  depend  strongly  on  the  blading  unsteady  aero- 
dynamic response. 

III.  RECOMMENDATIONS  FOR  FURTHER  WORK 

The  modeling  technique  developed  in  this  work  can 
be  used  beneficially  in  turbine  engine  compressor  design  and 
development  and  to  assist  in  solving  operational  problems 
which  may  arise  in  field  use.  There  are  both  improvements 
and  an  extension  which  can  potentially  increase  its  value 
further. 

Improvements 

Stage  characteristics . Stage  characteristic 
information  is  at  the  heart  of  this  modeling  method. 
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Efforts  should  be  made  to  obtain  improved  stage  character- 
istic information.  More  attention  should  be  given  to 
acquiring  accurate  stage  characteristic  data  during  com- 
pressor tests.  Theoretical  methods  for  computing  (or  at 
least  extending  the  experimental  results)  should  be 
developed . 

Unsteady  cascade  aerodynamics.  Continued  work  is 
necessary  in  this  area  to  define  the  unsteady  aerodynamics 
of  cascades,  particularly  near  and  in  the  stalled  operating 
regime . 


Crossflows . The  approximate  methods  used  for 
crossflow  in  the  model  are  somewhat  crude.  Specific  atten- 
tion should  be  given  to  improving  the  crossflow  represen- 
tation. A possible  approach  is  the  determination  of  a 
different  solution  method  to  solve  the  governing  equations, 
including  the  radial  and  circumferential  momentum  equations 
The  method  would  have  to  be  reasonably  economic  in  terms  of 
computer  storage  and  time  requirements.  This  would  also 
require  determination  of  a proper  representation  for  the 
radial  and  circumferential  forces  present  in  the  compressor 

Extension 


Post-stall  behavior.  The  model  could  be  extended 
to  compute  the  behavior  of  the  compressor  after  an  insta- 
bility was  entered  by  incorporation  of  unsteady  aerodynamic 
cascade  behavior  in  the  deeply  stalled  regime.  Because  of 
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the  highly  inefficient  processes  present  in  the  compressor 
in  deep  stall,  consideration  should  also  be  given  to  the 
metal-to-gas  unsteady  heat  transfer  processes.  Such  an 
extension  would  be  valuable  for  determining  methods  for 
returning  the  compressor  from  an  unstable  to  a stable 
condition. 
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a.  Influence  of  steady,  screen  generated  distortion 
on  turbojet  compressor  stability  limit  [4] 


b.  Influence  of  uniform  time  varying  flow  at  the 
inlet  of  a turbofan  engine  [5] 

Figure  4.  Typical  influences  of  external  disturbances 
on  compressor  stability. 
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c.  Influence  of  time-variant  distortion  on 

turbojet  engine  compressor  stability  limit  [6] 


Corrected  Airflow,  Ibm/sec 
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d.  Influence  of  time-variant  distortion  on  the 

low-pressure  compressor  of  a turbofan  engine  [7,  8] 

Figure  4.  (Continued). 
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e.  Influence  of  unsteady  and  nonuniform  flow  from 
the  low-pressure  compressor  on  high-pressure 
compressor  stability  limit  in  a turbofan  engine  [8] 


Corrected  Airflow,  -g— , Ibm/sec 

f.  Influence  of  nonuniform  total  temperature  on 
a turbojet  engine  compressor  [9] 

Figure  4 . (Continued) • 
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Figure  5. 


Surge  analysis  based  on  compressor  and  throttle 
curve  shapes . 
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Distance  Into  Conpressor 

Figure  6.  Representative  axial  flow  compressor  arrangement 
and  operation. 
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a.  Physical  system  modeled 


b.  Overall  system  control  volume 


c.  Elemental  control  volume 
Figure  7.  Physical  system  modeled  and  control  volumes. 
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a.  Rotating  fluid  element 


Figure  8.  Rotor  fluid  element  velocities  and  forces. 


152 


StaiiC  Equivalent 

Total-Pressure  Stage  Equivalent  Stage  Adlabatii: 

Ratio  Temperature  Ratio  Ef f Icioncy , n 


AEDC-TR-77-80 


1.2 

1 0 

0.8 

0.6 

0.4 

0.2 
1.28 

1.24 

1.20 

1.16 

1.12 

1.08 
1 . 5 

1.4 

1.3 

1.2 

C 

Figure 


6 0.7  0.8 

t,  $ 

itage  characteristics 


[46]  . 


53 


Figure  11. 


Oscillating  Foils, 
Steady  Stream 


Unsteady  airfoil  and  cascade  aerodynamics. 


AEDC-TR-77-80 


Ln 

OS 


k 


Figure  12.  Comparison  of  theoretical  and  experimental  unsteady 
cascade  lift  coefficients. 
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Figure  13.  Variation  of  dynamic  lift  coefficient  ratio  with  mean 
angle  of  attack. 
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Figure  14.  Unsteady  lift  coefficient  magnitude  ratio. 
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b.  Computational  planes 

Figure  15.  Computational  plane  designations  for  finite 
difference  solutions. 
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Figure  16.  Solution  procedure  for  one-dimensional,  time-dependent  model, 
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Figure  18.  XC-1  compressor  model  loaded  to  stall, 
N/Zff  = 90%. 
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a.  Imposed  nozzle  area  reduction. 

Figure  21.  NACA-8  compressor  model  loaded  to  stall,  N//0 
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b.  Compressor  pressure  ratio 


Figure  21.  (Continued). 
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Figure  22.  Comparison  of  model  computed  speed  character- 
istics and  stability  limits  with  experimental 
measurements,  NACA-8  compressor. 
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a.  Stages  1 through  4 

Figure  23.  Comparison  of  model  computed  flow  coeffici 
at  compressor  stability  limit  with  experin 
measurements,  NACA-8  compressor. 
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b.  Stages  5 through  8 
Figure  23.  (Continued). 
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Figure  24 . Comparison  of  model  computed  speed  character- 
istics and  stability  limits  with  experimental 
measurements,  J85-13  compressor. 
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a.  Inlet  total  pressure 


Figure  25.  Oscillating  inflow  loading  to  stall,  NACA-8  compressor, 
N//0  = 70%. 
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c.  Compressor  entry  airflow 


Figure  25.  (Continued) 
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g.  Flow  coefficients,  stages  7 and  8 


Figure  25.  (Continued) 
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Figure  26.  Stability  limits  for  planar  oscillating  inflow,  NACA-8  compressor. 
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a . Amplitude  ratio 

Figure  2B.  Variation  of  amplitude  ratio  and  phase  angle 

through  the  J85-13  compressor  with  oscillating 
compressor  entry  pressure  at  40  Hz. 
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b.  Phase  angle 
Figure  28.  (Continued). 
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Figure  29.  Dynamic  response  to  oscillating  compressor  entry  pressure, 
J85-13  compressor. 
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Figure  29 . (Continued) . 
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Figure  29.  (Continued). 
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a.  Amplitude  ratio 

Figure  30.  Variation  of  amplitude  ratio  and  phase  angle 

through  the  J85-13  compressor  with  oscillating 
engine  fuel  flow  at  40  Hz. 
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31.  Dynamic  response  to  oscillating  fuel  flow, 
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Figure  32.  Rapid  inlet  temperature  ramp  to  stall,'  NACA-8  compressor 
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Figure  32.  (Continued) 
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Figure  33.  Compressor  instability  caused  by  rapid  inlet 
temperature  ramps. 
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Figure  34.  Division  of  compressor  into  control  volumes. 
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Figure  35.  Control  volume  velocities  and  fluxes. 
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Figure  38.  Radial  distortion  effects. 


201 


12 


stall 


AEDC-TR-77-80 


c 


£ 


Angle  of  Attack,  a 


u 

<1 


m - 40 


c.  Unsteady  stall  effects 
Figure  39.  (Continued). 
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Figure  44.  Comparison  of  computed  and  measured  circumfer- 
ential distortion  attenuation  for  combined 
pattern . 
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Figure  44.  (Continued). 
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Figure  46.  Radial  profiles  with  uniform  entry  flow. 
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b.  Stage  2 

Figure  47.  (Continued). 
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Figure  47.  (Continued). 
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Figure  50.  Loss  in  stability  limit  pressure  ratio  with 
time-variant  distortion 
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b.  Hub  radial  distortion 
Figure  51.  (Continued). 
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Figure  53.  Model  computed  influence  of  circumferential 
pressure  distortion  on  stability, 

XC-1  compressor. 
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Figure  54.  Comparison  of  circumferential  distortion  model 
computations  with  experimental  results. 
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Figure  55.  Circumferential  crossflow  magnitudes  and 
effects. 
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Figure  59.  Comparison  of  circumferential  temperature 

distortion  model  computations  with  experimental 
results . 
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APPENDIX  B 

COMPRESSOR  DETAILS  AND  STAGE  CHARACTERISTICS 
Allison  XC-1  Compressor 

Figure  B-1  (Appendix  C)^  shows  the  overall  com- 
pressor map  for  the  Allison  XC-1  lift  engine,  four-stage 
compressor.  Key  dimensions  and  design  performance  are  given 
in  Table  B-1  (Appendix  C) Stage  characteristics  used  in 

I 

the  model  are  given  in  Figs.  B-2a  through  B-2d.  The  mean  of 
the  outer  and  inner  radii  was  used  in  the  stage  flow  and 
pressure  coefficient  calculations.  All  the  information  was 
taken  from  [43].  Dimensions  and  data  not  tabulated  in  the 
report  were  obtained  from  scaling  drawings  and  figures. 

NACA-8  Compressor 

Figure  B-3  gives  the  overall  performance  of  the 
NACA  eight-stage  research  compressor.  Principal  dimensions 
and  design  performance  are  given  in  Table  B-2.  Stage 
characteristics  are  given  in  Figs.  B-4a  through  B-4h.  Mean 
radii  were  used  in  the  flow  coefficient,  equivalent  pressure, 
and  equivalent  temperature  ratio  computations.  All  infor- 
mation on  the  NACA-8  compressor  was  taken  from  [44,  45,  46]. 


1 

Appendix  C. 


All  figures  and  tables  for  Appendix  B appear  in 
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General  Electric  J85-13  Coniprossor 

Overall  performance  of  the  J85-13  compressor 
modeled  during  this  study  is  shown  in  Fig.  B-5.  The 
variable  guide  vane  and  compressor  bleed  schedules  used  in 
the  model  (and  for  which  the  compressor  map  of  Fig.  B-5 
applies)  are  given  in  Fig.  B-6.  Key  dimensions  and  design 
point  performance  is  given  in  Table  B-3.  Rotor  inlet  air 
angles  for  which  the  stage  characteristics  were  calculated 
are  also  given  in  Table  B-3.  Stage  characteristics  are 
shown  in  Figs.  B-7a  through  B-7h.  The  first-stage  flow  and 
pressure  coefficients  were  re-defined  for  this  work  in  terms 
of  a semi-empirical  modified  flow  and  pressure  coefficient 
to  account  for  variable  guide  vane  effects.  The  modified 
coefficients  are. 


and 


m 


cos  3 

1 

— - sin  5 

ft 


(B-1) 


1 - sin 

L 

p 

where  0 and  ij;  follow  the  conventional  definitions  given  in 
the  text.  The  value  of  3 used  is  the  guide  vane  angle  from 
Fig.  B-6. 

Information  on  the  J35-13  was  taken  from  [47j. 


(B-2) 


■'ij 
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APPENDIX  C 

FIGURES  AND  TABLES  FOR  APPENDIX  B 
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Pressure  Coef 


a . Stage  1 

Figure  B-2.  XC-1  compressor  stage  characteristics 
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b.  stage  2 

Figure  B-2.  (Continued). 
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Pressure  Coefficient,  i|)  Efficiency 
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c . stage  3 

Figure  B-2.  (Continued). 
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Flow  Coefficient,  <J> 
d.  Stage  4 

Figure  B-2.  (Continued). 
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a.  stage  1 

Figure  B-4.  NACA-8  compressor  stage  characteristics. 
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b.  stage  2 

Figure  B-4.  (Continued). 
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c.  stage  3 

Figure  B-4 . (Continued). 
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Figure  B-4 . (Continued). 
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g.  Stage  7 

Figure  B-4.  (Continued). 


250 


Slago  Equivalent  StaKe  Equivalent 

Pressure  Ratio  Temperature  Ratio 


AEDC-TR-77-80 
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Figure  B-4 . (Continued). 
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Figure  B-7 . J85-13  compres 
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Figure  B-7 . (Continued), 
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Figure  B-7 
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Figure  B-7 . (Continued). 
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Table  B-1 . Detroit  Diesel  Allison  XC-1  Four-Stage 
Lift  Engine  Compressor  Model  Details 


Geometry 

Radius 

Stage 

Area 

Length 

Outer 

Inner 

(Station) 

(ft'^) 

(ft) 

(ft) 

(ft) 

1 

1.136 

0.2165 

0.6991 

0.3565 

2 

0.8173 

0.1558 

0.6903 

0.4651 

3 

0.6610 

0.1558 

0.6903 

0.5159 

4 

0-5045 

0.1386 

0.6838 

0.5540 

5 

0.4026 

- 

0.6816 

0.5800 

Design  Point  Performance 
N = 18,030  rpm 

= 38.7  Ibm/sec 

PR  = 4.92  (maximum) 
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Table  B-2.  NACA  Eight-Stage  Research  Compressor  Details 


Geometry 

Radius 

Location 

Area 

Length 

Outer 

Inner 

Station 

Stage 

(ft^) 

(ft) 

(ft) 

(ft) 

1 

it 

2.182 

0.475 

0.8333 

0 

2 

* 

2.182 

0.475 

0 

3 

* 

2.182 

0.475 

0 

4 

* 

1.679 

0.475 

0.4000 

5 

1 

1.679 

0.475 

0.4000' 

6 

2 

1.311 

0.450 

0.5265 

7 

3 

1.173 

0.283 

0.5667 

8 

4 

0.9960 

0.292 

0.6100 

9 

5 

0.8352 

0.258 

0.6546 

10 

6 

0.6860 

0.242 

0.6900 

11 

7 

0.5681 

0.250 

0.7167 

12 

8 

0.4737 

0.242 

0.7373 

13 

** 

0.4737 

0.242 

1 

0.7373 

14 

** 

0.3129 

0.500 

1 

0.7712 

15 

** 

0.3129 

0.500 

0.7712 

16 

Variable 

• 

0.8333 

— 

* 

Inlet  duct. 


* * 

Burner . 


Design  Point  Performance 
N = 13,380  rpm 

=65  Ibm/sec 


PR  = 10.0  (maximum) 
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Table  B-3 . General  Electric  J85-13  Eight-Stage 
Compressor  Details 


Location 
Station  Stage 


Area 


(ft^) 


Geometry 


Length 

(ft) 


Radius 


Outer 

(ft) 


Inner 

(ft) 


Rotor  Inlet 
Air  Angle 
(deg) 


1 

* 

1.300 

0.1856 

0.6708 

0.1900 

0 

2 

* 

1.300 

0.1856 

0.6708 

0.1900 

0 

3 

* 

1.300 

0.1856 

0.6708 

0.1900 

0 

4 

* 

1.150 

0.1856 

0.6708 

0.2812 

0 

5 

1 

1.066 

0.1856 

0.6458 

0.3156 

Scheduled 

6 

2 

0.8395 

0.1364 

0.3935 

5.4 

7 

3 

0.6821 

0.1116 

0.4472 

12.1 

8 

4 

0 . 5630 

0.09767 

0.4878 

14.8 

9 

5 

0.4690 

0.08725 

0.5175 

19.4 

10 

6 

0.4055 

0.07892 

0.5367 

23.2 

11 

7 

0.3698 

0.07750 

0.5472 

22.2 

12 

8 

0.3580 

0.09883 

0.5506 

21.6 

13 

** 

0.3542 

0.1000 

0.5506 

0 

14 

** 

1.312 

0.1000 

0.2683 

0 

15 

** 

0.3542 

0.1000 

0.5506 

0 

16 

** 

Variable 

- 

0.6458 

- 

0 

Inlet  duct. 

** 

Burner . 

Design  Point  Performance 
N = 16,500  rpm 

=43  Ibm/sec 

PR  = 7.7  (maximum) 
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NOMENCLATURE 

A Area 

ACC  Projected  area  of  control  volume  circumferential 

surface  in  the  circumferential  direction 
ACR  Projected  area  of  control  volume  circumferential 

surface  in  the  radial  direction 
ARR  Projected  area  of  control  volvune  radially-facing 

surface  in  the  radial  direction 

$ 

ARZ  Projected  area  of  control  volume  radially-facing 

surface  in  the  axial  direction 
AZ  Control  volume  area  normal  to  axial  direction 

a Acoustic  velocity;  blade  spacing;  a constant 

C Blade  chord  length 

C^  Blade  drag  coefficient 

Blade  lift  coefficient 

CF(r)  Stage  characteristic  correction  factor  for 

radial  variation 

Cp  Specific  heat  at  constant  pressure 

c^  Specific  heat  at  constant  volume 

CXFC  Circumferential  crossflow  coefficient 

CXFR  Radial  crossflow  coefficient 

D Blade  drag  force 

DLR  DyneUtiic  lift  amplitude  ratioed  to  quasi-steady- 

state  lift 

EPR  Stage  equivalent  pressure  ratio,  Eq.  (47) 
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ETR 

e 

F 

FC 

FR 

FZ 

f 

g(. . .) 

H 

HC 

HR 

HZ 

IMP 

IMPR 

i 

iiti 

ip 

j 

jm 


Stage  equivalent  temperature  ratio,  £q.  (46) 
Internal  energy 

Force  of  compressor  blading  and  cases  acting  on 
fluid,  including  wall  pressure  area  force 
Circumferential  force  acting  on  fluid  in  control 
volume 

Radial  force  acting  on  fluid  in  control  volume 
Force  of  compressor  blading  and  cases  acting  on 
fluid  in  the  axial  direction 
Frequency 
A function  of  ... 

A function  of  ... 

Total  enthalpy  flux 

Total  enthalpy  flux  transported  across  control 
volume  circumferentially-facing  boundary 
Total  enthalpy  flux  transported  across  control 
volume  radially-facing  boundary 
Total  enthalpy  flux  transported  across  control 
volume  axial-facing  boundary 
Impulse  function,  Eqs.  (5)  and  (160) 

Ratio  of  stage  exit  to  stage  entry  impulse 
functions 

Axial  location  index 
Adjacent  aixial  location,  i-1 
Adjacent  axial  location,  i+1 
Radial  location  index 
Adjacent  radial  location,  j-1 
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jP 

k 

L 

M 

ME 

MFF 

MFFE 


N 

0(.  . .) 
P 

PR 

PS 

Q 

R 

Re 

REI 

r 

rir.s 

s 

T 

TR 

TS 

t 

U 


Adjacent  radial  location,  j+1 

Reduced  frequency  based  on  full  blade  chord 

length.  Fig.  11 

Total  control  volume  length;  blade  lift  force 
Mach  number 
Equivalent  Mach  number 
Mass  flow  function,  Eq.  (74) 

Equivalent  stage  downstream  mass  flow  function, 
Eq.  (73) 

t 

Compressor  rotor  speed 
Order  of  magnitude  of  ... 

Stagnation  (total)  pressure 
Stage  total  pressure  ratio 
Static  pressure 

Rate  of  heat  addition  to  control  volume 
Gas  constant 
Reynolds  number 
Reynolds  number  index 

Radius;  coordinate  in  the  radial  direction 

Root-nean-square  amplitude 

Coordinate  in  the  circumferential  direction 

Stagnation  (total)  temperature 

Stage  total  temperature  ratio 

Static  temperature 

Time 

Torque  on  a fluid  element 
Axial  velocity 
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V 

V 

vol 

W 

WB 

WC 

WR 

WS 

WZ 


w 


X 

Y 

y 

z 


Velocity  components  other  than  axial  (designated 
by  subscript) 

Radial  velocity  component  in  three-dimensional 

model  development 

Volume 

Mass  flow  rate;  airflow  rate 
Compressor  bleed  flow  rate 
Mass  flux  across  circumferentially-facing 
control  volume  boundary 

Mass  flux  across  radially-facing  control  volume 
boundary 

Stage  shaft  work  added  to  fluid  in  control 
volume 

Mass  flux  across  axially-facing  control  volume 
boundary 

Circumferential  velocity  component  in  three- 
dimensional  model  development 
Energy  function,  Eg.  (8) 

Time  derivative  of  a volume  integral,  Eg.  (150) 

A general  dependent  variable 
Axial  coordinate 


Greek  Symbols 

a Angle  of  attack 

8 Flow  direction  angle  relative  to  axial  direction 

X Ratio  of  specific  heats 

d A difference 
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AP  Pressure  fluctuation  amplitude 

6 Ratio  of  compressor  entry  total  pressure  to 

standard  day,  sea-level-static  pressure 
e A constant;  an  error  magnitude 

Stage  adiabatic  efficiency 

6 Angular  location;  ratio  of  compressor  inlet 

total  temperature  to  standard  day,  sea-level- 
static  temperature 

6~  Angular  extent  of  low  pressure  region  in 

f 

circumferential  distortion  pattern 
X Blade  chord-to-axial  direction  angle  or  stagger 

angle 

p Density 

P 

tj;  Stage  pressure  coefficient 

T 

i;;  Stage  temperature  coefficient  (stage  loading 

parameter) 

u;  Angular  rotor  speed;  circular  frequency 


Superscripts 

Average  or  effective  value 

P Pertaining  to  pressure 

T Pertaining  to  temperature 


Subscripts 


0,1,2, 

ad 

avg 

crit 


Location  or  station  designators 

Adiabatic 

Average  value 

Critical  value 
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dyn 

exit 

H 

i 

im 

ip 

ijk 

j 

jm 

jP 

k 

km 

kp 

L 

loc 

m 

max 

min 

n 


P 

radial  avg 

ref 

rel 


Dynamic  value 

Pertaining  to  control  volume  exit  plane 
High 

Axial  location  index 

Adjacent  axial  location,  i-1 

Adjacent  axial  location,  i+1 

Three-dimensional  location  index,  refers  to 

value  at  location  i,  j,  k 

Radial  location  index 

Adjacent  radial  location,  j-1 

Adjacent  radial  location,  j+1 

Circumferential  location  index 

Adjacent  circumferential  location,  k-1 

Adjacent  circumferential  location,  k+1 

Low 

Local  value 

Modified  form;  minus  one 
Maximum  value 
Minimum  value 

Last  location,  e.g.,  in,  last  axial  location 
in  control  volxune 
Plus  one 

Average  value  along  a radius,  average  at 
constant  6 
Reference  value 

Relative  value,  usually  relative  to  rotor 
blade  motion 
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ring  avg 


stall 

ss 

sw 

wh 

z 


Average  value  along  a circle,  average  at 

constant  r 

Value  at  stall 

Steady-state  value 

Swirl 

Wheel 

Axial  component 
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